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Ionospheric Variations Associated with 
Geomagnetic Disturbances* 


I. Variations at Moderate Latitudes and the Equatorial Zone, 
and the Current System for the S, Field 


By Sadami MATSUSHITA 
(Geophysical Institute, Kyoto University) 


Abstract 


Variations in the densities and heights of the ionospheric Es 
and F2 regions during magnetic storms are studied. At the 
moderate latitudes in szmmer solstice and at the equatorial zone, 
the disturbance variation of the Es in both local time (S,) and storm 
time (Dsi) is obtained by statistical examination, although the 
individual values of the Es data may not always be accurate. 

The amplitude of the S, variation of the Es at the moderate 
latitudes is about 0.5 Mc/s in frequency (15% for mean ionized density), 
and the phase of the variation seems to be opposite to that of the 
F2. 

At the equatorial zone, the S, variations of the Es and the F2 
regions and the Dsf variation for f,F2 are rather peculiar for the 
phase, compared with those at the moderate latitudes. The geo- 
magnetic S, variation at this zone is also a little different from 
that in the moderate latitudes. 

On the other hand, the current system for the S, field is 
calculated by the dynamo theory, taking into consideration the aniso- 
tropic electrical-conductivity. A production mechanism of the Es is 
supposed, and the possibility that the ionospheric S, variations are 
due to the effect of a drift by the earth’s magnetic field and the 
electric field of the S, current is discussed. 


1. Introduction 


The morphology of the F2 variations associated with magnetic disturbances 
has been studied by many researchers [1][2][3][4][5]. But these variations at the E 
region, which may be an important region for geomagnetic variations, are not yet 
obtained. For the Sporadic-E (Es) region, intense occurrence of the Es during 
geomagnetic storms or bay disturbances at the high latitudes is reported by a few 
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researchers [6][7][8]. At the moderate latitudes and the equatorial zone, however, 
results of the disturbance variations of the Es have been obtained by the present 
writer [9], from statistical examination, although the individual values of the Es data 
may not always be accurate. Further statistical studies are made in this report, and 
some peculiarities of the disturbance variations for the Es and F2 at the equatorial 
zone, and those of the geomagnetic variation, are described. 

For the purpose of explaning the ionospheric disturbance variations in local 
time, the S, current-system must be considered. The idealized. current-system for 


(a) (b) 
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the S, field was investigated from geomagnetic data by S. Chapman [10], and more 
detailed examinations were carried out by M. Hasegawa [11], E.H. Vestine [12], L 
Harang [13], T. Nagata [14], and N. Fukushima [15]. The polar current of the S, flows 
parallel with the meridian of local noon in Chapman’s idealized current-system, 
but of about 09 h. in Hasegawa’s and Vestine’s systems, and of about 07.5h. in Nagata’s 
and Harang’s, as will be seen in Fig. 1. In the current system of the geomagnetic bay, 
the polar current flows parallel with the meridian of about 07.5h. according to H. 
Hatakeyama or H.C. Silsbee and E.H. Vestine, and of about 05h. according to N. 
Fukushima. On the other hand, the current focuses closed at moderate latitudes may 
be situated on the meridian of 06h. and 18h. for the average Sp variation, though 
it is a little different for the bay disturbance. 

As a theoretical approach to the S, current-system, T. Rikitake [16] has tried 
to explain the current system by the dynamo theory making the assumption that the 
electrical conductivity at the auroral zone is higher than that in the other regions. 
A further investigation of the dynamo theory for the S, current was carried out by 
N. Fukushima [17]. But both these cases were calculated postulating the isotropic 
electrical-conductivity, and cosequently the results of the calculation may not always 
be complete. 

Theoretical studies of the ionospheric conductivity have been made by D.F. 
Martyn [18], T.G. Cowling [19], M. Hirono [20] and K. Maeda [21]. The dynamo 
theory, taking into consideration the anisotropic conductivity of the ionosphere investi- 
gated by them, seems to succeed in interpreting the S, field. It has been reported by 
the present writer [9] that the Sp field is similarly calculated by the dynamo theory 
taking into consideration the anisotropic conductivity, under the assumption that the 
electrical conductivity at the auroral zone (and the polar cap) and at the equatorial 
zone is higher than that in the moderate region. A further investigation is carried 
out in this report. 

It has recently been reported by D.F. Martyn [5] and K. Maeda [22] that the 
Sp variation of the F2 layer is due to the effect of a drift of the earth’s magnetic 
field and the electric field of the S, current-system. But both studies for the Sp 
current-system are not entirely adequate. In this report, a theoretical basis for my 
statistical results using Martyn and Maeda’s principle is sought following a further 
examination of the S, current-system. 


2. The Disturbance Daily Variations (S,) 
2.1. The S, Variations of the Sporadic-E Region 


In the case of studying the variation of the Fs, the data may, unfortunately, 
not always be accurate, because the /Es value is influenced by incoming interference 
[23]. However, they are not without value. Accordingly, statistical examination must 


be conducted as far as it is possible. 
The S, variations of the Es have been studied statistically as follows: The 
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list of principal magnetic storms tabulated in the Journal of Geophysical Research and 
the report of magnetic storms issued from the Magnetic Observatory (Kakioka) have been 
prepared. Severe and moderately severe storms which were comparatively isolated, 
were selected from the lists and classified under three seasonal headings, Dec. sols- 
tice (Nov., Dec., Jan. and Feb.), equinox (Mar., April, Sept. and Oct.) and June solstice 
(May, June, July and Aug.). In order to derive the S, variation from the data of hourly 
values, the fEs value for the hour nearest to that of the beginning of the selected 
magnetic storm was picked up, and the monthly median value subtracted from it 
(4fEs). The same process was followed for the succeeding 48 hrs. (ie. lst and 2nd 
24hrs.) and for the preceding 24hrs. This process was repeated for each selected 
storm, and the result for each 24 hrs. was arranged and averaged in each local time. 
This method was applied to each season. 

When the resulting 4f/Es for the selected storms with sudden commencement 
is arranged and averaged in the 
storm time, the storm-time variation 
can be obtained, as mentioned in Para- 
graph 3. 

The average disturbance variations 
of 4fEs in local time, from these treat- 
ments of data at Kokubunji (139°29.3’E, 
35°42.4‘N) for 35 storms in the summer 
solstice of 1949-1952, are shown in Fig. 


Ist 24 hrs. 


2. In the same way, the variation for 
2004fEs/median fEs, showing the den- 
Argent sity variation in percentage, is also 

ay om plotted by the dotted line in the figure. 
— It will be seen in the figure that the 
Ist and the 2nd 24hrs. variations after 
the beginning of the magnetic storm 
have a greater diurnal variation than 
the preceding (-Ist) 24hrs. On the 
other hand, the 4/Es variation in local 
time during the ionospheric F2 storm’ 
0 6 12 18 24 (ie. the case in which the data are 


local time (hours) symbolized by the mark K) at Kokubunji 
Fig. 2. Average disturbance variations of AfEs : h : : 
in local time at Kokubunji for 35 storms mn, the, game period, 1s "shown at “mie 


in the summer solstice of 1949-1952, for bottom of the figure. The amplitude is 
preceding 24 hrs. (top) and succeeding 0.3M P 
48hrs. (middle two) from the beginning S c/s (10.7%) and the maximum 


of the magnetic storm. Bottom is AfEs occurs at 8.7h 9 
average-variation during the ionospheric é Tt re mie 
F2 storm at Kokubunji in the same median f. Es). , 

period. ‘The dotted line in percentage iati i 
Ee edit (Bs ' The S,(/Es) variation at Brisbane 


(153.0°E, 27.5°S) for 17 storms in the 
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Dec. solstice of 1949-1952 was also 
examined in the same way (see Fig. 3). 
The phase and the amplitude for all 
these results are shown in Table 1, 
and the 
storms are significant with a signifi- 
It will 
be seen in the table that the amplitude 
of the variation is about 0.5 Mc/s (15% 
for mean 


these amplitudes during 


cance-level of less than 109%. 


ionized density), and the 
phase is opposite to the S, variation 
of the F2 in summer. 

Reliable S,(/Es) variation at the 
moderate latitudes for the winter 
months and equinoxes was difficult to 
obtain, because of few occurrences of 
the Es in these seasons. 


As the S,(fEs) variation for the 
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BrisBANE 


-Ist 24 hrs. 


Ist 24 hrs, 


0 18 24 


12 
local time (hours) 
Fig. 3. Average disturbance variations of AfEs 
in local time at Brisbane for 17 storms 
in the summer solstice of 1949-1952. 


Table 1. 
Freee: 24hrs. before | Ist 24 hrs. after 2nd 24 hrs. after 
Stati Periods of t 5 Data the beginning of|the beginning of|the beginning of 
ations | data studied Scoendiie analized the storms the storms the storms 
study } ~ Fes ehh He lob ret aioe age 

Kokubunji| June solstice, Jo AfEs 0.14 Mc/s |16.5h.) 0.65 Mc/s |14.0 h.| 0.50Mc/s [08.0 h. 
1949-1952 

” 9. o 200 AfEs 4.6% |19.8h.| 19.0% 13.9 h.| 19.6% 06.9 h. 

median fEs 

Brisbane | Dec. solstice, 17 AfEs 0.32Mc/s |21.0h.| 0.76Mc/s |09.3 h.| 0.42 Mc/s |13.3 h. 
1949-1952 

_|Huancayo| Jan. 1951- 28 AfEs 0.38Mc/s |01.3 h.| 1.00Mc/s [23.3 h.| 0.46Mc/s |01.2 h. 


July 1952 


equatorial zone, the variation at Huancayo (75.3°W, 12.0°S) for 28 storms through all 
‘seasons from Jan. 1951 to July 1952 is shown in Fig. 4, and the amplitude (it has 
a significance-level of less than 10% during the storm) and the phase will be seen in 
Table 1. The variation maximum in the equatorial zone seems to occur at earlier 
hour than that in the moderate latitudes. 

Although the individual values of the h’Es data may have less accuracy and 
the period examined was short, the S,(h’Es) variation at Brisbane was also 
‘studied in the same way, for 7 storms in the Dec. solstice, 1951-1953. The results 
are shown in Fig. 5, and the amplitude and the phase are compared in Table 2. The 
former during the storm have a significance-level of less than 15%, but the variation 
must be further examined for longer period data. 
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Km. 


6 -lst 24 hrs. 


-Ist 24 hrs. 


SSCS Pm LS 6 ist 24 hrs, 
Ist 24 hrs. 
0 
6 
2nd 24 hrs. 
Huancayo 0 
0 6 fe 18 24 
local time (hours) Fis 
F - SB 
Fig. 4. Average disturbance variations of AfEs 6 SBANE 
in local time at Huancayo for 28 storms 
in Jan. 1951-July 1952. 0 6 12 18 24 


local time (hours) 


Fig. 5. Average disturbance variations of Ah’Es 
in local time at Brisbane for 7 storms 
in the summer solstice 1951-1953. 


Table 2. Sp(h'Es) at Brisbane 


i 


24hrs. before the beginning of | 1st 24 hrs. after the beginning | 2nd 24 hrs. after the beginning 
the storms of the storms of the storms 
P, / ty | P, | ty | Py : ty 
1.9 Km. | 16.5 h. | 2.2 Km. | 01.2 h. | 3.3 Km. | 23.9 h. 


The S,(h'Es) variation for Huancayo was also investigated during the same 
period as the study of fEs. P, and t, for the 2nd 24 hrs. are 1.7Km. and 56h. 
respectively, although they seem to be less accurate, tested statistically. 


2. 2. The S, Variations of the F2 Layer at Huancayo 


The Sp variations in the critical frequencies, S,(f,F 2), and in the minimum 
virtual heights, S,(h'F2), of the F2 layer at Huancayo are studied for the same 


period as the treatment of fEs, following the same method as in Paragraph 2.1. The 


S)(//F2) variation is obvious, as will be seen in Fig. 6, though the S,( foF2) is not so 
clear. The amplitude and the phase of the variations are shown in Table 3, and the 


former for the S,(h’F2) has a significance-level of less than 10%. The variation 


Ionospheric Variations Associated with Geomagnetic Disturbances LL5 


maximum seems to occur in the day 
time for the S,(f)F2) and at mid- 
night for the S,(h'F2). This mode is 
also somewhat peculiar compared with 
that of the moderate latitudes [1][3][4] 
[5], as in the case of the Es. 

For the purpose of a further 
Ist 24 hrs. study, the S,(h'F2) variation at Huan- 
cayo was computed, by subtracting, 
for each month of 1940-1944, the 
average daily variation for the five 


‘-Ist 24 hrs. 


international magnetically quiet days 
from the daily variations on the five 
To international disturbed days. The 
average disturbance variation of h’F2 
in local time during 1940-1944, obtained 


from this treatment, is shown in Fig. 


Huancayo 
- 
7.a, and the harmonic results for each 


0 6 12 18 "4 year and for each season are shown by 


local time (hours) the harmonic dial in Fig. 7.b. Taking the 
Fig. 6. Average disturbance variations of Ah’F2 


in local time at Huancayo for 28 storms mean value of these harmonic results, 


in Jan. 1951-July 1952. the amplitude is 14.2Km. and the 
Table 3. S,(F2) at Huancayo 


24 hrs. before the 1st 24 hrs. after the 2nd 24 hrs. after the 


iiiet doit sgn Sth af the beuinning, of the yee = the 
ae i ee a P| ty 
Af F2 0.29 Mc/s 12.1h. 0.35 Mc/s 14.2h. 0.27 Mc/s 15.4h. 
Ah’F2 3.0Km. 18.9 h. 20.0 Km. 00.8h. 20.0 Km. 01.3h. 
23 J 1 


Huancayo 


local time (hours) 


Fig. 7.a. Average disturbance variation of Ah'F2 20Km. 


0 10 
c ; i = ig. 7.b. Harmonic dial of the Sp(A’/F2) varia- 
in local time during 1940-1944 at Huan- Fig tiptataaacaoneri tachi vee ot 
i 1940-1944 and for each season (Dec. 
solstice, equinox and June solstice) 
during the period, 
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maximum occurs at 01.7h. in local time. These values resemble the preceding results, 


but reliable S,(foF'2) variation could not be obtained from these treatments. 


2. 3. The Geomagnetic S, Variations at Huancayo 


As mentioned in Paragraph 1, the current focuses of the average Sp current- 
system enclosed in the moderate latitudes, seem to be situated on the meridian of 
06h. and 18h. in local time. Since the S, variations of the ionosphere at Huancayo 
differ from those in the moderate latitudes, the S, variation of the horizontal compo- 
nent of geomagnetic variation was also examined for 5 years, 1940-1944, by the method 


mentioned in Paragraph 2. 2. 


‘ MEAN 
Huancayo Fig. 8.b. Harmonic dial of the Sp(H) variation 


Pte tr eer gas at Huancayo for each year of 1940- 


0 6 2 18 24 1944 and for each season (Dec. sols- 
local time (hours) tice, equinox and June solstice) during 


; é ae the period. 
Fig. 8.a. Average disturbance variations of 


geomagnetic horizontal component in 
local time at Huancayo for five years, 
1940-1944. 

The results for each year and the mean variation of the five years will be seen 
in Fig. 8.a, and the harmonic dial of Fig. 8.b shows the harmonic results for them and 
for each season. Taking the mean in the dial, the maximum occurs at 02.4h. in local 
time, with an amplitude of 17;. It follows that the center of the current in the 
equatorial zone shifts to a three or four hours earlier meridian than that in the 
moderate latitudes, and that the current density is a little larger there. ; 

It may be concluded, accordingly, that the Sp variations at the equatorial zone, 
for the Es region, the F2 layer and geomagnetic variation, have some peculiarity com- 
pared with those in the moderate latitudes. 


wor" 
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3. The Storm-Time Variations (Dst) 


The method employed in deriving the Ds¢ variations has been described in 
Paragraph 2.1. Following that method, the data at Kokubunji, Brisbane and Huancayo 
were examined for magnetic storms with sudden commencement during the same 
period as the study of the S, variations. As the variation was somewhat irregular, 


mean values for every six hours were obtained (dividing a day into four periods) to 
see the mode of variation. 


BRISBANE 


KoxkuBUNII 


BrisBANE 


Huancayo 


Huancayo 0 1 2 
storm time (days) 
Fig. 10. The Dst(h’Es) variations for Brisbane 
(4 storms in the Dec. solstice of 1951- 
1953) and Huancayo (19 storms in Jan. 
"1951-July 1952). 


1 Huancayo 

7 storm time (days) 

Fig. 9. The Dst(fEs) variations for Kokubunji 
(20 storms in the June solstice 1949-1952), 
Brisbane (9 storms in the Dec. solstice 
1949-1952) and Huancayo (19 storms in 
Jan. 1951-July 1952). 


From these studies for the Dst 
(fEs) variations at three stations, fEs 
seems to increase during several hours 
after the sudden commencement of 
magnetic storms, as will be seen in 


i 0 1 2 
lk ‘ storm time (days) 
The Dsi(h'Es) variations were Fig. 11." The Dst(f°F2) and Dst(h'F2) variations 
also examined for Brisbane and Huan- at Huancayo for 19 storms in Jan. 1951- 
July 1952. 


cayo. A small increase of about 2 Km. 
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on the 1st day seems to occur in both moderate and equatorial regions (see Fig. 10). 

Fig. 11 shows the Dst(f)F2) and Dst(h'F2) variations for Huancayo, and both f,F2 
and h’F2 seem to increase during magnetic storms. This result of the increase of fyF2 
is the reverse of that in the moderate latitudes in summer, and it coincides with the 
result obtained by Appleton [1], though the method used is not the same. 


4. The Current System for the S,-Field 
4. 1. The Aurorai Zone and the Equatorial Zone of the Es Region 


For the purpose of calculating the current system for the S, field, the 
position and breadth of the auroral and the equatorial zones must be decided at the 
outset. 

It has been reported by the writer [24] that the world-wide distribution of the 
fEs may be divided into three groups: the auroral zone (intense at night), the 
equatorial zone (abnormally intense in the day time) and 
the intermediate zone (proportional to cos x, where 7% is 
the sun’s zenith distance). From a further study for mon- 
thly median values of fEs (CRPL-F series report) at the 
auroral and the equatorial regions, the distributions in Fig. 
12 and 13 were obtained. 

Fig. 12 shows fEs distributions for three seasons on 
geomagnetic latitudes at night (mean values from 22 h. to 
02h.). It increases at latitudes higher than about geo- 
magnetic 60°N, and reaches the maximum at about 
69°N. In still higher latitudes, the fEs decreases, and at 
Resolute Bay (geomagnetic 83°N) in Canada the Es becomes 
very small. The auroral zone of the Es region, accordingly, 
may be considered as a region from 65°N to 75°N in 
the geomagnetic latitude. 

As will be seen in Fig. 13, fEs in the day time 
(mean values from 09h. to 15 h.) is abnormally large near 
the magnetic equator. This abnormality resembles the 
range distribution of daily variation of the horizontal 
magnetic force by J. Egedal [25] in the left part of the 
figure, and also resembles the range distribution observed 

40 60 80 °n at South America, India and Togoland by J. Egedal and 
Fig. 12. oop aa iat i A.A. Giesecke [26]. Accordingly, this intense Es zone may 

three seasons (Dec. be due to the effect of the Hall drift [20] related to the 

solstice, equinox and intense eastward current flow that produces the large daily 


June solstice, 1951) 
on geomagnetic magnetic variation, as mentioned in Paragraph 5. As 
latitudes at night. another peculi ‘ f 

(Tier ob axe t-vagus peculiar character of the Es at this zone, jEs 
in the day time.) values at Huancayo may have some correlation with the 


Sun-spot cycle (the correlation coefficient between them 
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being 0.60 for monthly median values and 0.66 for overlapped means), as will be seen in 
Fig. 14, although such a relation can not be obtained in the moderate latitudes. Judging 
from these results, the equatorial zone may be considered as a region between 5° and 


—5° of the magnetic dip. 


4. 2. The Dynamo Theory, Taking into Consideration the Anisotropic 
Electrical-Conductivity 


The electrical conductivity of the ionosphere and-——taking this into considera- 
tion-—-the dynamo theory have recently been studied by many researchers [18][19][20] 
[21]. These results will be briefly introduced below. 

Ions and electrons are distinguished respectively by the suffixes 7 and e. The mass, 
charge (in e.m.u.) and number density of 7th molecules are denoted by ,, e, and 7,, 
respectively, and the charge of an electron is specially denoted by —e. Thus, the 
formulae for the ionospheric conductivities are expressed as follows: 


the direct conductivity is do= 7D} (n,/M,v,), 

the Pedersen conductivity is isi ey] (n,/m,){v,/(v,?+o,2)}, 

the Hall conductivity is 02 = (n,e/ A){o?/(v.2+,")—,7/(v2+0;7)}, 
and 030; +02"/01, 


where, w,=¢,H/m,, H is the intensity of the geomagnetic field, and v, is the colli- 
sional frequency of the rth molecule with the neutral one. 


KM It must be noticed that o, is larger than o, at 
ae 105+25Km. (see Fig. 15) and that o, has a pro- 
= nounced sharp maximum at about 100 Km. level near 
180 the magnetic equator. 


When x- and y-axes are taken respectively as 
southward and eastward directions, the effective 
direct conductivities are 


Gxx> 0y0;/Q), 


Fyy = (G90, Sin?g + o105 cos*¢)/Q, 


ie 1 ea 10 aU as and the effective Hall conductivities are 


Fig. 15. Variations with height = i 
Oxy = —Oyz = 0902SiN i 
of o1e/Me (1), o1/m% (2) wate Pee tae 
and o:/n, (3). where ¢ is the magnetic dip angle, and 


(after M. Hirono) Q =o sin’?¢+oe cos*¢ 
— f: . 


If the permanent magnetism of the earth is assumed to be a dipole magnet 
with its axis coincident with the earth’s rotation, then, the upward component of H, 
W348 
fi=Cens 0, 
2 
where C= — 3 (€.m.u.) and @ is co-latitude. 


Now, we shall consider a case where the region between /, the bottom of the E region 
and the upper height h,, makes a uniform horizontal tidal oscillation with the velocity 


et Ng 
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potential 
; ¢ = knPm (cos @) sin (mA+an), (1) 


where Pr(cos@) denotes associated Legendre function of Neumann, and 2 is longitude 
reckoned from the midnight meridian, and k® is taken to be positive. When E" is 
the dynamo field (=C, H, and C, is the mass velocity), and EF is the electrostatic 
field having three components X (southward), Y (eastward) and Z (upward), then the 
electric field EF is their sum. 

First we shall consider all the regions except the zone near the equator. In the 
dynamotheoretical calculation, the effect of the region higher than h#, must also be 
taken into consideration. This time we can use the same horizontal electrostatic 
field for both the lower (Ap<h<h,) and the upper (%;<h<h,.) region of the iono- 
sphere, as a first approximation. The conductivity }} and current density J integrated 
with height for these upper and lower regions are 


hy hg 


pi=|orda Sty" =[opda, (b=, Wy, 4, YX 1,2, 3) 


ho hy 
P=\Jdz, 9 =\JTdz. 


Then we get the following equations, 
[’= Sii2(0H, +X) +217,(—uH,+ Y), 
L,- = >¥4.(0H,+ X) + 314,(—uH,+ Y), 
Fr = 3 eX Hy Y, 
dy” = 2. -A+ 359° Y; 


for the lower region, 


for the upper region, 


where R, is the radius of the conducting sphere, and 7=0//R,00, v=0¢/R,sin 904. 
The components of J=/’+J” are accordingly, 


I, = 32,-0H, 4+ (ee te) X — by UH + (Qi + Dey) Y, 
y= 34,-0H,4+ (242+ 309) X Py 5, 4H, 4+ (ey 4 D7) Y. 


(2) 


We adopt 2,=180 Km., /,=300 Km. and take the gas density distribution assumed 
by M. Hirono [20] from rocket observations. Then, except near the equator, 
pias =)}},/sin*¢, Dw = Dh, Dey =)),/sin dg, for both regions, 


and 
$32=686X10"%; Spex eh Le *5 


Se 159x10"° << , 2. 2h” =0. 
Therefore, we can neglect the terms containing }Jp” in Eq. (2), and substitute 1 for 
the denominators sing and sin’, as the first approximation. Then we get from Eq. 


(2), for the northern hemisphere, 
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f= D(H, + X)+>3(- uH, + xe 


L, as —)>).(vH, +X)+3)i(— uH, Swe Y). 
(Suffix Z for >} is neglected.) 


(3) 


When we introduce the current function R, I,=0R/R,sin9-01, I,= —OR/R,00. 
From these equations and 0X/04—0(sin9-Y)/00=0 and also Eq. (3), X and Y are 


eliminated, giving 


OR 0 OR 


‘ a J 
sin Joe tart oa (sin? 3) = Rokk) 


O(vH.) 4 “ar en) | 
OA 00 7 


where K=}};+(d}2)?/3. This equation is satisfied in the range 0<0<6;, and, from 
the statistical result mentioned in Paragraph 4.1., 6,=85° is adopted. The solution 
of this equation is 


m m mn ™. 6 $ m ™. Y A m | 
R= GahomED Clik { Pitas Costan | sin (md +oft)+Cnatan™--c05 mi +o) |, 


where Cn; and Cn are constants for each value of m. In the same way, we get the 
solution for the zone 85°<0<z. 

In the zone 85°<0<95°, we can assume H.~0. Thus it follows that o,,,=0, 
Fyy~O3, d,y=0. We take o,, as infinite because of d)>03. Then X=0 or ss = 0, and 
we have 


i,=K,«¥,. «where .. K~3ig (4) 


Accordingly we get 0R/00 = K,(0S/sin 004). Electrostatic potential S is eliminated, 


giving or (See) =0. On the assumption that K, is independent of 9, the solu- 


tion of this equation is 


ae ¥ 0 m = 
R= (m+1) @m-+1) Chm K log (tan z) {Cnssin (mA + Gn) + CmsCos (mad +an)} ? 


where Cn; and Cn, are also constants. 
The current function R, for the S, field is accordingly, 


; ° ° a m iy m™ 0 
in the zone O°<0<85°, R= pm 4 Cli: K| { Pitart Coatan x 
Sin (mA + a7") -+ Cg tan™$-cos (mA +- af) | ; 


in the zone 85°<0<95°, R,= 


m ns 0 
(m+1) @m+l) mK log (tan 5) (5) 
{Cm sin (A +0") + Cm cos (mA+am)}, 


in the zone 95°<0<180°, = im) Ona FD CRA Ki [{ Pliss—Cmtan> | 


sin (mi-+ ait) —Cnstan™ 9 cos (m’-+ 0 |. 


Electrostatic components X and Y for each zone can be obtained from Eqs. 
(3), (4) and (5). 


eS ee 
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For m=1, 
in the zone 0°<6<85°, 


Ck;' cosé 1 1 
gr oy Soler 2 eet GRAAL RAD EME 
i | 


Ee 6 Thest (TGQ, — ta) (a +r) |, 


“ae hs es | 
! jee & ttcost (C+rC)|sinQ@+as')+) 5 cos@ 
i 1 


+§ TpeoaT C.-C) jeos(a-ta |, 
in the zone 85°<0@<95°, 


1 Ch K 
re 6 Re K, cosec# {Cssin (A +ay!) + Cscos(a +a}, 


in the zone 95°<6<180°, they are obtained in the same way, where y=)>}:/>)1. 
For m=2, 


in the zone 0°<6<85°, 


Ck,” 
x= R. eG —sin26+Fe 5 (pasty Cit 1G) c0s 2a-+a5% 


15 (1+cos 0)? 
2 6 ; 
+forsinaacnre af Bp lre cine 0 
Chee? 4 sin 
t cloy ig R, Fe {esind— 15 ressoy (Cit1C,)}sin +a") 


0 
+ {2rsin20-+75 apeer TCs cy cos(2i-+ax")], 


in the zone 85°<0<95°, 
A= 
eo 


z 
oS os cosec 0} Cysin (22-+0s")+Cycos (22+ a"). 


The constants C;, C., C; and C,, in the current functions and the electrostatic 
components, are determined from the continuity of the normal component of the 
current and the tangential component of the electric field along the boundary. They 
are as follows, assuming K,=2K, 


for m=1, 
C, =(0.0226—0.00067?)/1"1,. C2=0.00737/T"1, 
(6) 
3 = —0.5134/11, C,= —0.07657/T1, 
where J°,=0.0929+0.00217?, 
for m=2, 
C, =(1.0726—0.09677)/P'2, C2= 0.65107 /T"s, 
(7) 
C3= —24.152/T2, C,= —6.254r/Ts, 


where J’,=0.9321+-0.06257°. 
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4. 3. The S, Current-System by the Dynamo Theory, Taking into Consideration 


the Anisotropic Electrical-Conductivity 


The S, current-system may be calculated as follows by the dynamo theory, 
assuming that the electrical conductivity at the auroral zone (and the polar cap) is 
higher than that in the moderate latitudes. 

As will be seen in Fig. 16, the electrically conductive layer of radius FR) is 
divided into seven parts, which are founded on the statistical results obtained in 

Paragraph 4. 1. The height-integrated 
I a 0" 15° effective-conductivity in each region is 
assumed as follows: 
zone I, northern polar cap, 
0°<@< 15°, conductivity AK, 
zone II, northern auroral zone, 
15°<0< 25°, conductivity BK, 
OWT 85° zone III, northern middle pes 
95° 25° <@< 85°, conductivity K, 
zone IV, equatorial zone, vee’ 
85° <@< 95°, conductivity 2K, 
zone V, southern middle zone, 
95°<0<155°, conductivity K, 
VI a re - 165° zone VI, southern auroral zone, 
VI ~—-180° 155°<0<165°, conductivity BK, 
Fig. 16. Division of zones I, II, III, IV, V, VI zone VII, southern polar cap, 
sai 165° <8<180°, conductivity AK, 


where @ is colatitude, A and B are constants larger than unity and K= Dh +(2)?/Dh.- 

At about 100 Km. level, , is much larger than’ o,, as already mentioned in the 
preceding paragraph (see Fig. 15). Accordingly, in consideration of the increase of 
the conductivity at the auroral zone (and the polar cap), it may be assumed that only 
DJ: increases, and >}, does not change. Then, the height-integrated Hall conductivity 
at the polar cap and the auroral zone may be A’>}, and B’>}, respectively, where A’ 
and B’ are constant. In the preceding paragraph, we put PAP />) as y, and now the 
ratio at the auroral zone becomes B’y, and the ratio at the polar cap A’y. It may be. 


estimated, accordingly, that the height-integrated effective- -conductivity at the auroral 
zone has the relation 


BK=B{S+(3))'/S3 7 
=>) {1+(B’r)"}. (similarly for A and A’) 


Now, we can obtain the equations of the current functions and electrostatic 
fields at each zone, for the velocity potential ¢ of Eq. (1), following the same principle 
as in Paragraph 4. 2. The current functions are 
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in zone I, 


FE m 6 : 
R= oT Om tI) CAB { Pits +Cmstan™ 5 lsin md +0 


— Costan”-5-cos (ma a om | ’ 


in zone II, 
= eo ve bd ™m. iy a4 
®t Demt1) Chi: BK] | Pitsa+Costan sy + Cnscot™ 5 z)* in (md +o) 


+| +? Cns tan” = + Son te cos (ma + wa] ; 


in zone III, 


m 


a 7] pe 
R= (m-+1)(m-+1) K ie mei Cutan" + Cnscot”-- > sin (mA + af) 


fe { Cns tan” as 4. Cais cot” 5 } cos (mi + o| 9 


in zone IV, 


m 


n= (m+1)(2m+1) mK log (tan 5 x) {Cmy; Sin (md +- 0%) + Cry. cos (mA +an)} « 


In zones V, VI and VII, R is expressed in the same way as in zones III, II and 
I respectively. Twelve constants in these equations are determined by the 
conditions that the normal component of the current and the tangential component 
of the electric field are continuous along the boundaries. 

As a first approximation, we consider that the current function is divided into 
two parts, namely R=R,+Rp,, where the current system derived from Rp corresponds 
to the S, field, and that derived from R, causes the S, field. The equation for FR, is 
already discussed in the preceding paragraph. Accordingly, the remaining part Rp is, 
in zone I, ' 

1 me ; m 
Ro= GED emTD eR ((A-D Pirsit (ACms— qm) tan > }sim (onl ae 
0 
+ (ACmng = qm) tan”--cos (ma + a) | “ 


in zone II, 


1 ea “x! 1 ikke 
Rp = (m-+iyam+1) k| {(B —1)Pis1 + (BCnz —qm,)tan >) + BCn,cot 2 sin(ma +om 


Fa {@ Cog —gm2) tan™-5 + BCne cot” 5 leos (ma +02) , 


in zone III, 


a ; mo “>| i A m 
Ro= Gap nemEt HA | {(Cor—amadtan” 3 + Cmscot™ x jsin (mel + ot) 


+ { (Cng—Qme) tan" + Cmocot-5-}c0s (mA +0) | ' 


in zone IV, 


a ah Soe) Pi A i Cmig—Qgms)cos(mi+amn)}; 
Ro= peep em ry OE oR (tan 5 x) (Conis—qns) Sin (mA +084) + (Cmsa—Gns)COS( ) 
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where q,, is C,, in the equations of R,. In zones V, VI and VII, R, is expressed in 
the same way as in zones III, II and I respectively. The components of the electrostatic 
field X, and Y, for the S, in each zone are also obtained in the same way, as shown 


in the succeeding paragraphs. 
4. 4. Results of Numerical Calculations for the Diurnal Wind Variation (m=1) 
As a simple case, m in Eq. (1) is assumed to be 1, namely 
=k, P,'(cos 9) sin (A +4,'). 


Then, the current function Rp and electrostatic components Xp and Y, in each zone 
are, 


in zone I, 


1 
Ry=GCh'K| ((A-DP:'+(AC—q)tan-p-}sin(@ +0;')+(AC,—q,)tan- cos (2 +a:)], 


sae ale ee cos(A+a;') + 15 (A’— 100820 
+ + AG Ge-patel snares] 
Yo=g| FA fi htsin (ht at) +{ (A'—1)cos9 
4. + ages Foss md eek 
in zone JJ, 


eysing +a,') 


+ {(BC.—q:)tan7-+BCycot-7 loos +a], 


Ry=+-Ch, 'K| {(B- I)P2'+(BCy~q,) tan $-+-BC,cot-$- 


a RG (rea OG Joos tas) 
es (B'—1)c0s20+-4 (FES Spare Bret ee) sin@+a'], 
Yo TR | 8 (Ped et eo ina) 
+HF@ —1) cosd+-h B (en ot eg pons tat | 


1 
Rp a howal fc qi)tan-s-+ Cecot 5 \sin(i +a!) 


+ (Coa) tang + Cio cot | cos(t-+a19) i 


__Ch;! 1(G+rco— Yea Pg TCiw—C, 


OF ir ® ee Co—rqitge _ 7Cs+C. ° 
IG ( 1+cos@ T—cos8 ) sin(i-+- a], 
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_Ck'f 1 (CrtrGQ—q—re , rCwo—C.\ . 
Yo= eel I+cos9 Tt tad) sin +a) 


1 ( rCr—Co—rqitqe 


Cs + Cro 
6 1+cos 0 


‘3 1—cos0 costa) |, 


-E 


in zone IV, 
1 7] : 
Ry= =| Chi'Klog (tan +) {(Cy1—gs) sin (A +ay")+(Ci2— qs) cos (A+-a, yes 
xX,=—0 


1 Ck, : 
Yp= — = R, Fe e086 8 {(Cyx— 4s) 8in (A405) + (Crx—qs) 608 (2-04), 


where q:, gz, gs and q; are given as C,, in Eq. (6). In zones V, VI and VII, they are 
obtained in the same way as in zones III, II and I. 

Numerical calculations for many sorts of values of y, A’ and B’ are carried 
out, and a typical example showing the mode of those results is given below. 

When we adopt >}, =6.86 x 10-°, }},=1.45x 10-8, as mentioned in Paragraph 4. 2., 
then y=2.114. At that time, the calculated current-systems for the following four 
examples, 

1! M=f, “hence, A=1, B’=2 hence B= 3.45, 
2, A‘=15 hence A=2.02, B’=2- hence. B=_3.45, 
oA =sd hence (Asi, B'=6 hence B=29.60, 
4, A'=3 heace. A=754, B'=6. hence. B=29.60, 
are given in Figs. 17. a, b, c and d, where a,' is assumed to be 114°, and the density 


II 


of stream lines in Figs. a and b is expressed as twice that in Figs. 17. c and d. 
From these examples, the following conclusions may be obtained: 

1) When the height-integrated Hall conductivity >}, at the auroral zone 
increases to only twice that of the other zones, as will be seen in Fig. 17. a, the 
focuses of the current circulating through the middle and the equatorial zones are 

situated on the meridian of about 04h. 16h., that is about 4 hours later compared with 


(a) (b) 


CH 
aii 


a 


ZA 
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(d) 


<_) 


—— 


Fig. 17. Calculated Sp current-systems for m=1 and y=2.114 (viewed from the north 10% 
a) A’=1 and B’=2, b) A’=1.5 and B’=2, c) A’=1 and B’=6, d) A’'=3 and B’=6. 


the calculated S, current-system in the case of isotropic conductivity [17]. With 
further increase of >}, at the auroral zone (or B’), the focuses shift more anticlock- 
wise till about 06h. 18h. meridian. It may be found by comparing Fig. 17. a with 
Fig. 17. ¢: 

2) The direction of the polar current hardly changes for any values of 
B’, as shown in both cases of Figs. 17. a and c, but it turns anticlockwise with 
the increase of >}, at the polar cap (or A’), in so far as B’ is constant. Fig. 17. 
d is the case where the increase of >}, at the polar cap is a half of that at the 
auroral zone (A’/B’=0.5). If A’ increases almost to B’, the polar current then finally 
flows as in Fig. 17. b,—it does not become the return current of the auroral current. 

3) The electric current of 3.9210? Amp. in Figs. 17. a and b, (and its double 
in Figs. 17. c and d) flows between the successive stream lines in the direction 
indicated by arrows, assuming that the current intensity for the S, current-system 


Table 4. 
Case ‘} 2 3 4 Chapman’s 
Pou eount of tN | 9.85 |) 3.35 hee ee 27.50 
ps cane pt ae 1.11 0.532.891.9958 22.50 
gangs seem 2.44 2.82 «| 4.27 4.26 5.00 


(unit in 104 Amp.) 


agrees with Chapman’s analysis. The amount of current is compared in Table 
4, for the total auroral zone current, the polar closing current (the current 
circulating through the auroral zone and the polar cap) and the middle zone closing 
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current (the current circulating through the auroral, the middle and the equatorial 
zones). The current intensity of these results is rather small, compared with that of 
Chapman’s Sp, current-system. But excessive increase of B’ for the purpose of 
producing a larger current may be impracticable. 

4) The direction of these calculated currents does not coincide with the Sp 
current estimated from the geomagnetic variation as shown in Figs. 1. a, b, c, and d, 
for any values of y, A’ and B’, but it seems to be in reverse phase. If a,! can be 
assumed to be 294° instead of 114°, the difference in phase of both current systems 
dissapears, as may be supposed in Fig. 17.d. Accordingly, if the wind blows reverse 
direction to the wind producing the S, current, the S, may be explained by this 
principle, though the current intensity rather small compared with Chapman's. As 
an hypothesis, it may be reasonable, for our purpose, to consider double layers [27] 
[28] wherein the air motion is assumed to be reciprocally inverse and the conductivity 

. differs, or to assume some coexistence of such a wind system during the storm with 
the wind producing the S, current. 

5) At the equatorial zone, the phase of the calculated S, current is about two 
hours earlier than that in the middle zone for both the cases a,1=114° and a,!=294°, 
as the 0-line in each figure shows. This result seems to coincide with the preceding 
statistical result for geomagnetic horizontal component at Huancayo. 


4. 5. Results of Numerical Calculations for the Semi-Diurnal 
Wind Variation (m=2) 


Assuming that m=2, then ¢=h,"P.*(cos@)sin(2A+o,"). Rp, Xp and Yp are 
also calculated in the same way. They are, 


in zone I, 
: 6 : 
Ro= jg Che | [(A-DP+ (AC a) tan’ |sin 22+ a") +(AC.—q.)tan*-7-cos(22-+a2") |, 
_ Chef 4 2 pW GtATG.—q—-14: 2 
Ap= SE yen ee cos (2A+ az ) 
/ pers me F 
+{24/—1)r sind @eos'0—1) +35 ae ie oe et | sin@a-+a.) |, 
Chl 4 eg GtArG—a—1a.: p 
Yo=-Fe | —y5sind (1-+c0s0)? sin (2A +-a,”) 
/ = so 
+ [2¢Aa’ apr sin20 + 75 sind ST Pte cos earrac’) |, 


in zone II, 


2 Y] 2 Y 7 2 
Ro=jgCheK| ((B-DP#+ (BC;—q:)tan’-> + B Cycot? 5-}sin(22 +a") 


Ch2f 4. C,+-B'rC: =f tgs _ Bir Co ) : 
ey Ry [asin¢(“S* areas (i—cos 6) ) °° (24 +42") 


’ he ByCs—Cs—raitage  BirCs+Ce\ .; 2 
+ {2@e’—rsino@eos0—1)+755in0 ( (iF c0s6)! “— T—cos 8), sin (24+ a”) |, 
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ee FET — 35 8i00( Cr+ BI Cs— Qi —142 Bree) sin (24+ 2”) 


Is 15 (1+cos 6)’ (1—cos 0)? 


: 4. pf BrOs—Cs—rmtag: , BrOstGs 2 
+ {2(B'—rsin20 +35 sind ( (L-Fc0s8)? 24 (1—cos 8)? cos (24+ az”) |, 


in zone III, 


6 a es 
Rp = 5 Ck!K it (Cr— qu) tan’ +Cz cot" } sin (24 +a") 


+ {cs — qe) tan? y+ Cro cot? >| cos(2i-+a) |, 


x, [ ey (Se (oa) cos (2A +42”) 


FURS Ib © (1-+cosé (1—cosd)? 

4. 1Cr1—Co—71 Ut 42 _ pet os) in (2h »| 
+75 sind ( (1+cos6)? (1—cos 6) sin(2A+as') |, 

Che a x Cr+rCo—qi—1qe2 ew) : 2 

Yo= Fe | — igen? (“i pooety + G—eoagy inl tax*) 

be 1C,—Cy—rqitqe pot x.) 2 
+7gsind (J “(i-cosby * (i—cosdys) Cosa ae, 
in zone IV, 


Rp =— Ch,°K log (tan >) {(Ciu—qgs) Sin (24 +02") +(Ci2—qu) cos(24 +42*)}, 


Xp=0, 


ahs Ce + cosec 9 {(Cy,— qs) sin (24 +42") +(Ci2—gs) cos (24+ a,")} , 
0 e 


ey yates 


where qi, g2, qs and q, are given as C,, in Eq. (7). In zones V, VI and VII, they are 
expressed in the same way. 
If we assume in this case that >}, =4.8x 107°, >},=1.44x10-8, then y=3. At that 
time, the calculated current-system for the following two examples, 
1, A’=1 hence A=1, B’=6 hence B=32.50, 
2, A’=3 hence A=8.20, B’=6 hence B=32.50, 
are given in Figs. 18. a and b, where a," is assumed to be 292°. 
The density of stream lines for the polar closing current in Fig. 18. a is 
(a) (b) 


Fig. 18. Calculated Sp current-system for m=2 and y=3 (viewed from the north pole) 
a) A’=1 and B’=6, b) A’=3 and B'=6, 
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expressed as half of that for the middle zone closing current of the figure and for 
Fig. 18. b. If the current density for the S, current in this case is assumed to be 
62,000 Amp., it may be estimated that the electric current of 5.46x10° Amp. (and 
double that amount, for the polar closing current in Fig. 18. a) flows between the 
successive stream lines in the direction indicated by arrows. The amounts of the polar 
closing current and the middle zone closing current for Fig. 18. a are 9.83 x 10! Amp. and 
2.50 x 10? Amp. respectively. 

The second example seems to be too complicated at the polar cap, but in the 
first example of Fig. 18. a, it may be noticed that the amount of the polar closing 
current is about four times that of the middle zone closing current, and the ratio in 
Chapman’s Sp, current is 4.5. 

The direction of the current in Fig. 18. a is also worthy of notice. As men- 
tioned in Paragraph 4. 4, if the wind blows in the opposing direction to the wind 
producing the S, current, the Sp, may be explained to a certain degree. But when 
the wind is the same for both the S, and the S, currents, this case of semi-diurnal 
wind variation may possibly explain the S, current, from the point of view of the 
direction of the current. In such a case, the daily variation of electrical conductivity 
must be taken into consideration. 

Dynamo theoretical calculation for these conditions can be treated by means of 
a perturbation method [29]. It may be possible to assume that the height-integrated 
effective conductivity in each region is 


in zone I and VII, K=AK, ( The conductivity does not change. 
in zone II and VI, K=B&K,(1+¢cosA), The conductivity at midnight is 
(1+¢)/(1—¢) times that at noon. 


The conductivity is proportional to cos x, 
in zone III and V, K=K,(1—¢’sin9cos4), [| in the case when the sun’s north decli- 
K=2K,(1—¢’sin 9 cos2) nation is 0, where x is zenith distance 

=2K,(1— : 


i me IV, 
RO of the sun. 


where ¢ and ¢’ are smaller than 1. The current function R is assumed to be 

R=R,4+Rp, 
in the same way as in the preceding paragraph. 
Then, Rg= Rot O Roy , 

Rp=Roa to Roy tER' vy ’ 

where ¢R’>, is a perturbation term corresponding to BK,¢cosi, and both ¢’R’,,, and 
¢/R’p., are perturbation terms for —Ky¢’sin9cosi. Ryo and Rp, are the expressions 
for R, and R, when ¢=¢/=0. Thus, we can calculate the Sp current for the semi- 
diurnal wind variation with a daily variation of electrical conductivity. Though these 
results will be discussed in the second report, it will be easily seen that the 
current at night in Fig. 18. a predominates, and accordingly it becomes similar to the 
Sp current obtained from the geomagnetic variation, 
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5. Ionospheric Variations due to the Drift 


Ionospheric variations due to the effect of a drift by the electric field and the 
earth’s magnetic field were discussed by many researchers, and it has been suggested 
recently by D.F. Martyn [5] and K. Maeda [22] that the S,(F2) variation is also due 
to a drift effect. Following their principles, some theoretical approach to explain the 
statistical results is carried out below. 

For positive ions in the E layer, the equation 

On, /0t=Q—an,n,—0(n,W,)/02z, 
will hold, where Q is ion production per c.c. per sec., a is the effective recombination 
coefficient, and W is vertical upward velocity. The terms for horizontal velocity are 
neglected, since they are estimated to be small. When A( = *§ — 1) is assumed to be 


independent of height, we have 
0n,/0t=q—an—0(n,W,,)/02, (8) 


where g=Q/(A+1). 
According to M. Hirono [30], W, is expressed as 


W, =F (K,-cos$-sing-J.+K, -cos $-Jy), 


where 
K,=—D{—M?-L,1,4+M(L+L,)}, Ke=DthitLh+M(LL,—D}, 
D=L,/{(L2+*L,)1+M?*)}, M=L,(L,+«L;)/(L2+«L,), 
L1= | o:/¥ |; L,= | ,./». |; 
cK=2A+1. 


Although the second term of the right hand side of the equation is smaller than the 
‘first in all zones except the equatorial, J, is very large near the E level in 
‘the equatorial zone when 4 <1 [31], corresponding to the intense geomagnetic daily 
variation at this zone. It follows that W, is much larger there than in any 
other zone. Accordingly, in consideration of Eq. (8), some abnormality in the electron 
‘density or the height of the ionosphere at about the E level in the equatorial zone 
may be supposed. In the E layer, it seems to be difficult to obtain statistically such, 
an abnormality. In the Es region, however, fEs is abnormally large in the equatorial 
zone as shown in Paragraph 4. 1., and the breadth of the zone coincides with that of 
intense J, zone, statistically and theoretically. 

Another expression of W by Hirono is 


* RAVES 3 
bis W,= 6r(S, — Pre (Ee +-X)+(de-+ Per Ev +¥)s (9) 
. Pos (l /m,v,){cos’¢+v,2/(v,2-+0,2)+sin’?¢d}, S,=(1/m,»,)sind-cos¢-a,2 /(V,2+ Wy2)>, 
dr =(1/H) cos $-0,?/(v,?+0,), P=e'>}\n,P, =0,Cos*¢ +o, sin’¢, 
T= M0 (4. —4) = 7020S g, S= ey} NS, aa (a —d;) sin od *COS d. 


The variations of the E layer due to a vertical drift associated with the Sy cho the 
L fields can be discussed by Egs. (8) and (9). The S)(£) variation may be discussed 
in the same way for the S, field, taking X, and Yj, instead of E,'+X and E,'+Y 
in Eq. (9). At lower altitudes than 120Km., the first term of the right hand side of 


Ionospheric Variations Associated with Geomagnetic Disturbances 133 


Eq. (9) is much larger than the second term in latitudes higher than about 5°, but it 
is very small near the magnetic equator. For the phase of the S,(E) variation, 
therefore, Y> is effective in the equatorial zone and X, is important in other zones. 

In the E layer, however, it seems to be very difficult to obtain Statistically the 
Sp variation and lunar tide, as in the case of the lack of abnormality at the equatorial 
zone. This is because a is about 1078 cm’/sec., so that the second term in Eq. (8) decreases 
the drift effect. In the Es region, on the other hand, intense lunar tidal variations 
were obtained by the writer [82], and the S, variation is also obtained as already 
mentioned. Moreover, fEs is abnormally large in the equatorial zone. Therefore, 
these results suggest that the Es is much affected by the vertical drift, and in 
other words, they suggest a formation mechanism of the Es. That is to say, it may be 
supposed that the Es is a comparatively thin Jayer formed of ionic clouds, having 
the electron density of a few or several times that of the normal E, and that the 
intense ionic clouds are due to a drift effect on patches or blobs [33] of ionization, 
of different density from their surroundings, caused by meteors or other things. At 
the auroral zone, however, the Es may be formed mainly by the charged particles 
emitted from the sun. 

Thus, it may be estimated that the S,(Es) variation may be due to a drift 
effect associated with the S, field, and that Y, may be effective for the drift in the 
equatorial zone and X, may be important in other zones, for the purpose of con- 
sidering the phase of the S,(fEs) variation. In the case of ,=2.114, A’=3, B’=6 
and a,!=294°, the minimum of X, in the middle zone occurs at 07-08h. in local time, 
and maximum of Y, in the equatorial zone occurs at 05h. Accordingly, already 
mentioned statistical results that the S,(fEs) has a ‘maximum in the forenoon at the 
moderate latitudes, and that there is a phase difference between the middle zone 
and the equatorial, can be explained by these considerations. 

In the F2 layer, ionospheric variations due to the effect of the vertical drift 
may be obvious, as a is small. For this layer, W can be expressed by the second 
term of the right hand side of Eq. (9). In the case of the S,(F2) variations, therefore, 
Y, may be important in considering the phase of the variation, and Y, has a phase 
difference between the middle and the equatorial zones. The increase at about 05h. 
in the S,(h’F2) variation for Huancayo (see Fig. 6 and 7.a) may be explained, when 
Y, at this zone is considered, as mentioned above. A.A. Weiss [34] has pointed out 
that the greater the vertical velocity the smaller the electron density of the F2 layer 
in the day time, and K. Maeda [21] has also obtained the same conclusion, Accordingly, 
it may be supposed that the S,(f,F2) is at a minimum in the forenoon at the middle 
zone and that there is some phase difference between the equatorial and the middle 
zones, as obtained statistically. 

It may be concluded, therefore, that the vertical drift associated with the Sp 
current differs in the middle and the equatorial zones, and its effect is opposite for 
the F2 layer and the Es region,—that is to say, the density decreases in the former and 
increases in the latter. These matters will be examined further in the second report. 
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6. Concluding Remarks 


1) The Sp and Dsé variations in the densities and heights of the Es and F2 
regions are studied. The amplitude of the Sp variation of the Es in summer solstice 
at the moderate latitudes is about 0.5Mc/s in frequency (15% for mean ionized den- 
sity) and the phase is opposite to that of the F2. 

2) At the equatorial zone, the ionospheric S, and Dst variations are 
a little peculiar compared with those at the moderate latitudes. The geomagnetic 
S, variation at this zone is also a little different—that is to say, the amplitude seems 
to be larger and the maximum occurs about three hours earlier than that in the 
moderate region. 

3) The current system for the S, field is calculated by the dynamo theory, 
taking into consideration the anisotropic electrical-conductivity, on the assumption that 
the conductivity at the auroral zone is higher than that in the other zones. In the case 
of the diurnal wind variation, it may be possible to consider double layers wherein 
the air motion in one layer is assumed to be the inverse of the wind producing the 
S, current in the other and the conductivity differs between the two, or to assume 
the coexistence of such a reverse wind system during the storm with the wind 
producing the S, current. In the case of the semi-diurnal wind variation, the daily 
variation of the electrical conductivity must be taken into consideration. The latter 
case seems to give rather better results, judging by the current density for the S, 
current-system estimated from the geomagnetic variation. In both cases, the peculiarity 
of the geomagnetic S, variation at the equatorial zone can be explained by these 
calculated current-systems. 

4) The ionospheric Sp variations may be due to the effect of a vertical drift 
by the earth’s magnetic field and the electric field of the S, current. From the 
equations of vertical velocity and an estimation that the production of the Es is due 
to a drift effect, the statistical results that the phase of the S, variation is opposite 
between the #2 and the Es, and that there is a phase difference between the 
equatorial and the middle zones in both ionospheric regions, may be explained. 
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Latitude Effect of the Intensity Decrease in 
the Cosmic-Ray Storm* 


By Sekiko YOSHIDA and Yoshiko KAMIYA 
' (Physical Institute, Nagoya University) 


Abstract 


The latitude effect of the cosmic-ray decrease associated with 
the magnetic storm was determined from the data obtained in 13 
periods of observation at various stations over the world. The result 
was in accordance with Nagashima’s electric field theory. 


Nagashima‘? proposed an explanation of the world-wide decrease (47) of the 
cosmic-ray intensity associated with the magnetic storm, assuming an electric field, 
and calculated the latitude effect of 47 according to his theory. The latitude effect 
was found by Forbush“ and by Wada“ through the analysis of observations. But, 
neither the number of stations nor the procedure in their analysis is satisfactory to 
compare with the theoretical result. We have made this analysis, taking the follow- 
ing points into account. 

1) In order to decrease the error due to the variations except the cosmic-ray 
storm, only the storms with large decrease are subjected to.the analysis. The crite- 
rion adopted by us is 4/,>2%, where 4I, is the difference between the averaged 
intensity during the three days just before the storm and the daily mean intensity on 
the day of intensity minimum, calculated from the data observed at Cheltenham. 


Among the materials ready to our hands, 13 storms were available by this criterion, 
as shown in Table I and IL. 


Table I Cosmic-ray storms 


N Period (U.T.) used in Date of S.C. (U.T.) Ah 

oO. , ‘ A 
the analysis of magnetic storm % 

1 Apr. 16—May. 10, 1937 Apr. 24, 25, 26 3.3 
2 | Jan. 10—Jan. 31,°1938 Jan. 16, 17,* 22, 25 5.8 
3 Apr. 20—Apr. 28, 1939 Apr. 23, 24 2.9 
4 May. 3—May. 11, 1939 May. 5 3.0 
5 Aug. 9—Aug. 20, 1939 Augy) 11, 12)), 2.4 
6 Mar. 27—Apr. 5, 1940 Mar. 29 2.0 
7 Feb. 26—Mar. 9, 1941 Mar. 1 2.6 
8 Feb. 26—Mar. 6, 1942 Mar. “1 6.0 
9 Feb. 1—Feb. 14, 1946 Feb. 7 6.1 
10 May. 1—May. 20, 1946 May. 6 2.0 
11 Jul. 26—Jul. 30, 1946 Jul. 26 6.3 
12 Jan. 24—Jan. 29, 1949 Jan. 24 3.0 
13 Apr. 12—Apr. 15, 1949 Apr. 12 3.0 

* Gradual commencement 

ii sn Nee e 


* Contribution of Geophysical Papers dedicated to Prof. M, Hasegawa on his sixtieth birthday. 
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Table II Cosmic-ray stations 


A f : 

meee Station Longitude PR Altitude | Equipment Reference 
H Huancayo 75.3°W 0.6°S 3350m |I.C. model C meter] (5) (7) 
Ec East China Sea| 124 E 22 N sea level |I.C. Neher type (6) 
Fy Tokyo 139.8 E 25.4. N sea level /1I.C. Neher type (7) 
Ty I.C. Steinke type (7) 
Ts I.C. Nishina type (6) (7) 
Ts C.T. +85° (7) 
Te Teoloyucan 99.2 W 29.7 N 2285 I.C. model C meter} (8) (9) 
Ca Canberra 1441 E 45: 3S 800 I.C. volume 17.67 (10) 
Ha, Hafelekar 11.3 °E 48.4 N 2300 I.C. Steinke type (11) 
Hag CT: (11) 
Ch Christchurch 172.6 E 48.6 S sea level |I.C. model C meter | (5) (8) (13) 
F Friedrichshafen| 9.3 E 49.0 N sea level |C.T. (7) 
C Cheltenham 76.8 W 50.1 N sea level |I.C. model C meter] (5) (7) (8) 
B Boston 2” Ww 54. ON I.C. model C meter| (14) 
M Manchester 0 56.5 N sea level |C.T. +60° | (7) (15) 
Be Bergen 5 E SL = ON sea level |C.T. (16) 
G Godhayn 53.5 W 79.9 N 


sea level I.C. model C meter] (5) (8) 


C.T.: Counter telescope ~ LC.: Ionization chamber 


2) In order to decrease the error due to the local slow variation, such as the 
atmospheric effect, especially the temperature effect, only a short period including 
“each storm was subjected to the analysis (see Table I). 

3) In order to decrease the error 
due to the variation in _cosmic-ray 
anisotropy during the storm,“ obser- 
vations at various. stations were 
arranged in the universal time (U.T.), 
and the simultaneous ?4-hour averag- 
ed intensities were treated as a set of 
data. 

As an example of the procedure, 
the day-to-day plot of a cosmic-ray — 
storm is shown in Fig. 1. The correla- 
tion diagrams between the decrease at 
a certain station and that at Chelten- 
ham are shown in Fig. 2. From each 


- of diagrams, the correlation coefficient 


? ee , "1946. 
ite ites feat alae July. sa vy, the constant factor k, and the pro- 
observed at various stations. 


bable error of k (small sample) were 
calculated. and shown in Table III together with the results from the other storms. 
The factor & may be considered as the ratio of the decrease (4/) at a certain station 


138 S. Yosmipa and T. Kamiya 


Table III (uncorrected for altitude effect) 


H Ec ly we 
7 PRR yl 7 t ee r 
1 |.96/0.99+.04 .97 
2 |.9010.78-+.07|.86|0.87-+.10)2|.94/0.73+.06 
3 |.94|0.84+.09 
4 |.97|0.73+.06 
5 |.93/0.81+-.08 
6 .83 0.88+.15 3 .72 0.80-+.10) 
7 |.7610.89+.17 
8 |.9710.93-+.02 3}.89]0.91+-.08].97 
9 |.99/0.80-+-.03 
10 |.83/1.08+.10 
11 |.98/1.08+.07 1}.99/0.88-+.05 
12 |.96/1.00+-.07 4|.97/0.95+-.06 
13 ‘|.98/0.79+.08 4|.99]0.82+..06 
mean|  |0.89+.07 0.85-4.05 


to that at Cheltenham. &-values or 
their average obtained at each geomag- 
netic latitude are shown in Fig. 3. The 
averaged point is the weighted mean 
of k-values obtained at the same latitude 
on individual storms, were the weight 
is the accuracy of individual k-value. 
The sign of error attached to the 
averaged point is the unbiased disper- 
sion among the individual k-values, the 
weight being taken into account. The 


curves are the result of Nagashima’s 


% 


theory,™ corresponding to his two 


Fig. 2 Correlation of cosmic-ray decreases 
observed at pairs of stations, in July, 
1946. that the latitude effect is in accordance - 


with Nagashima’s electric field theory, even if the result of high altitude stations 
(Huancayo, Teoloyucan, Canberra and Hafelekar) are omitted. There are sea level 
stations at about the same geomagnetic latitudes as those of high altitude stations 
respectively, except Huancayo. The altitude effect obtained from the comparison of 
these pairs of stations are listed in Table IV, together with the values expected from n 
Nagashima’s theory. The observed effect is not appreciably different from the theo- 
retical result. For the four high altitude stations, the correction for the altitude effect 
‘were done according to Nagashima’s theory, using the average of the two values 
obtained from his extreme assumptions respectively, and both of the corrected and 
uncorrected values are shown in Fig. 3. The corrected values and the sea level 


observations make a nearly smooth curve of the latitude effect. Thus, Nagashima’s 


extreme assumptions. 


From this figure, it can be seen 
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Cheltenham : k=1 


we Ca Ha Ch F B M Be G 
k r k ix?) Sak i ok mes Hi teak eee rile eR r k 
_ 1.35-£.07).87|1.13 +-.08|1|.98|1.30-4-.05].91|1.19-+.09 
2|.91/1.29-+.05 
.95/0.91+.04|1|.91/1.41-++.09|.99|0.57-+.02 .94|0.97-+.06 
.98]0.76-+.05 .94/0.96-+-.09 .97\0.94 +.09 
.97|0.76-+.06 .93]0.41+.09 84/1.15-+.16 
.85)0.91+4.13 .29|0.99+.67 .96|0.95+.06 
.92/1.48+.12 90/L.11+.11 .93)1.28-+.09 
9010.97 +.10 
1.10+.02 .90)0.67+.08),97/1.134-.03 .99}1.05-+.04 
.97|1.08-+.04 
89]1.67+.08 
92)1.154.15 99}1.06+.06 
-85)0.83-4.15 .78|1.02+.19].94|1.02+.21 

.94|0.96-+.15 

1.18+.03} |0.944.15| | |1.33.03| |0.924.21| | it | {1.06.07} ) L12+.13 


No. : see Table I 


Fig. 3 Latitude effect of the 
cosmic-ray decrease. 

@@O Average of k-values 
obtained at the same 
latitude 

e k-value obtained for a 
single storm 


O High altitude station, 
uncorrected 
@ Corrected.for the alti- 
tude effect 
| Normalization point 
(Cheltenham) 
Table IV Altitude effect 
Stations Altitude effect 
Geomagnetic high altitude sea level | observed theoretical _ 
latitude (h) (s) kn/Rs assumption I. assumption IT. 
30°-25° Teoloyucan (2285)m| Tokyo 1.39+.06 1.30 1.15 
48 -50 Hafelekar (2300) | Cheltenham 1.33+.03 1.45 1.30 
45 -50 Canberra ( 800) | Cheltenham 0.94+.15 LedS 1.09 


EARS EE —E———EE—————————EEE—EE nae 


electric field theory may be considered to be able to explain the latitude effect and 
also the altitude effect of the intensity decrease in the cosmic-ray storm. 

This work suported by the scientic research found of Ministry of Education. 
The authors wish to express their sincere thanks to Dr. J.W. Beagley of the Christ- 
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The World-Wide Variation of Cosmic Ray Intensity by 
the Electro-Magnetic Field* 


By Kazuo NAGASHIMA 


Earth Science, Kyoto University 


Abstract 


The variation of cosmic ray intensity by the static electric field 
surrounding the earth was calculated in a previous paper. In it a 
greater advantage was indicated in explaining the cosmic ray intensity 
decrease at the time of the magnetic storm by the electric field 
hypothesis than by the magnetic field hypothesis suggested by S. 
Chapman. In the present paper, the intensity variations of the 
ionizing and neutron components produced by the variation of the 
static electric field are calculated in detail extending over various 
altitudes and latitudes. By comparing these calculations with observa- 
tions, it can be concluded that this electric field hypothesis is possible : 
in explaining the phenomenon of the world-wide variations in intensity 
pointed out recently by various authors.“ 

It may be possible that the diurnal variation of cosmic ray 
intensity can be produced by an electric field not derived from 


potential. 


® 


Introduction 


Up to date, various types of cosmic ray intensity variation have been reported. 
Among these, the present paper is concerned with the following: 
[A] The world-wide variation of cosmic ray intensity (not dependent on local time). 
This involves the following phenomena: 
(a) The world-wide decrease of cosmic ray intensity at the time of the cosmic 
ray storm. 
(b) The world-wide fluctuations of cosmic ray intensity. 
[B] The diurnal variation of cosmic ray intensity. 
In order to explain phenomenon [A], any theory or hypothesis must correspond the 
following conditions : 
[I] The variation is world-wide. 
[II] The latitude dependence of the variation. 
Based on his experiments,“ H.V. Neher pointed out that: 
(1) The changes in the ionizing component at sea level are nearly independent 


of latitude. 


* Contributed to Geophysical Papers dedicated to Prof. M. Hasegawa on his sixtieth birthday. 
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(2) The changes in the ionizing component increase with the increase of 
latitude. 
[III] The altitude dependence of the variation. 

(1) The fluctuations in the ionizing component at high altitudes are larger than 

at intermediate altitudes or at ground level.“ 

(2) The changes in the neutron component are several times those of the ioni- 

zing component.@> (23) 

It was pointed out in a previous paper“ that the variation of electric potential 
on the earth causes such world-wide variations in cosmic ray intensity, and suggested 
that the comparison between the variations at different stations would be a key 
towards determining the true cause of the production of such a world-wide variation. 
Recently, according to this electric field hypothesis, W.H. Fonger explained the pheno- 
menon (2) in [III]. But, to explain all the above mentioned phenomena by this electric 
field hypothesis, it is necessary to calculate the cosmic ray intensity variations, pro- 
duced by the electric field, over various altitudes-and latitudes, more precisely. 

In § 1, the general formula of the cosmic ray intensity variation which is produced 
by the static electric and magnetic field is discussed. Generally, in such a static 
electric and magnetic field, the cosmic ray intensity variation is due not only to the 
variation of the primary differential intensity spectrum caused by the variation of the 
static electric field but also to the variation of the cut-off energy (see § 1) caused by 
the variation of the static electric and magnetic field. In § 2, neglecting the influence 
of the variation of cut-off energy, there is calculated the intensity variation which is 
produced by the variation of the static electric field. The variation of cut-off energy 
cannot be determined unless the state of the static electric and magnetic field sur- 
rounding the earth is known. In case of a simple model of the field, the variation 
can be obtained analytically and is discussed in App. I. But the altitude curve 
of the intensity variation caused by the variation of cut-off energy can be obtained 
without such a model of the field and is discussed in § 3. In § 4, comparisons 
of the intensity variations calculated theoretically in § 2 and § 3 are made with those 
obtained by observation. 

The diurnal variation of cosmic ray intensity is a world-wide phenomenon which 
is dependent on local time. The solar magnetic field theory was advocated by M.S. 
Vallarta” to explain this phenomenon, but at low latitudes this theory seems to be 
inadequate. In § 5, the possibility of explaining this phenomenon by the electric field 
hypothesis is discussed, although the origin of the electric field cannot be pointed out. 


§ 1. The General Formula of Cosmic Ray Intensity 
Variation by the Static Electric and Magnetic Field 


Let N(A, x) denote the directional cosmic ray intensity, observed at geomagnetic 


latitude 4° and atmospheric depth xgcm72, and the following relation holds at the 
normal state ; 
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NA, *)=3) M4, 2(E, #)t4,2(E)GB, -i0+.. Cl} 
AJE), 4,2 


where 7z4,,(E) is the primary differential intensity spectrum of the nucleus, its atomic 
number A and its charge number Z, and E is the kinetic energy per nucleon. 
my, 7(E, x) is the number of secondary particles at xgcm~* produced by one of the 
primary cosmic rays whose kinetic energy is E per nucleon and called the over-all 
multiplicity. F, 4°, denotes the minimum kinetic energy per nucleon of primary 
particle which can be injected upon the earth at 4° from a given direction and is called 
the cut-off energy. 

Strictly speaking, the formula of Eq. (1) is not correct because it may happen 
that the primary spectrum might in some way be interrupted in some energy region, 
and consequently N(A, x) cannot be expressed by only the upper and lower limits of 
integral in Eq. (1). For example, such a situation occurs when we intend to explain 
the diurnal variation of cosmic ray intensity by the solar magnetic field.“ In such 
a case, however, Eq. (1) can be easily expanded. The case chiefly discussed here is 
one where Eq. (1) holds. 

If the state of the static electric and magnetic field surrounding the earth varies 
from the normal state, then, by the influence of this variation, the primary cosmic 
ray, reaching the earth, loses its energy to an amount of 4E per nucleon, where 


— Z eee ee eeae 
4E=7- cA. (2) 


4@ is the variation of the earth’s electric potential relative to the normal state. In this 
case, Eq. (1) changes as follows, 


N(A, x, 4E, 6E,)= p>) i m4, AE, x[1—-L(E+ 4E, 4E))-t4, (E+4E)dE,:+++(3) 
4 YE) 4, 2492), AZ : 
where L(E, 4E) is Liouville’s effect®*“ and given by Eq. (4). 


2(E+me')4E—(4EY (4) 
(E+ me’? — (moe) ” 


LE, 4E)= 


where mc? is the rest energy per nucleon. In Eq. (3), 0F,4,z is the variation of 
cut-off energy, which is caused not only by the variation of the magnetic field but 
also by the variation of the electric field and cannot be uniquely determined unless 
the state of the electro-magnetic field is determined. dE, in the parenthesis on the 
left hand side of Eq. (3) denotes the occurrence of the variation of cut-off energy. On 
the contrary, the variation of the primary spectrum in Eq. (3) is determined if only 
the electric field is derived from potential“) and is independent of the form of the 
static electric and magnetic field. More generally, the static magnetic field which is 
not symmetrical to the earth’s magnetic dipole axis could change the primary spect- 
rum in the low energy region. For example, we see that the solar magnetic field 
produces such a situation if we consider that, as pointed out above, the solar magnetic 
field has an agency to interrupt the primary spectrum in some energy region, and 
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moreover that the primary particles exist in some degree in such an interrupted 

energy region of the primary spectrum because, as pointed out by H. Alfven“’? and 

and J.A. Wheeler et al.“ there exist the primary particles in the forbidden region 

of the solar magnetic field. The intensity variation produced by such a static magnetic 

field depends on local time, so we will discuss this in § 5, and not consider it here. 
The intensity variation relative to the normal state is defined by 


(A, x, 4E, 8E)= NA = SAE OO xe oa SS 


As seen from Eqs. (5) and (3), y is not proportional to 4E.0Y 

As mentioned above, the causes which produce the variation y, are divided into 
two classes as follows: (1) the variation of primary spectrum which is produced by 
the variation of electric potential and (2) the variation of cut-off energy which is 
produced by the variation of static electric and magnetic field. Hereafter, an examina- 
tion is to be made of the variation y by these two causes individually. 


§ 2. The Cosmic Ray Intensity-V ariation by the Static Electric 
Field (not considering the Variation of Cut-Off Energy) 


Recently, the existence of heavy particles in the primary cosmic rays has been 
pointed out.“ But, as their intensity spectrums in the high energy region are not 
yet well known, we assume that all the primary particles consist of protons only and 
adopt as its intensity spectrum that produced by H.V. Neher.@>(® 


" 0.048 
MEME + 00E a te (E=E,=0.8Ber,) | 


=0, (E <E£,) 
where E£, is the energy, in the lower energy region than which the primary particles 
do not exist, and is called hereafter the general cut-off energy to avoid being confused 
with the cut-off energy E,. This spectrum leads to a smaller absolute value of 
intensity compared to that observed by others.“ But, it only is necessary that the 
distribution of the primary spectrum is correct and it is not necessary that the 
absolute value of intensity is correct, because the variation y in Eq. (5) depends only 
upon the distribution, but not upon the absolute value. In this point, it is reasonable 
to adopt the above spectrum. 
; Neglecting the variation of cut-off energy, and considering only the vertical 
cosmic ray intensity, Eqs. (1) and (3) becomes as follows, 


N(A, x)= Vem (E, *)*2(E)dE, ++00050 (7) 
N(, x, 4E, 0)= Vem (E, )[I1-L(E+4E, 4E)]-i(E+4E)dE, «++++++- (8) 


where E, denotes the cut-off energy of cosmic rays which are incident upon the earth 
at 4° from the vertical direction. To calculate y in Eq. (5), it is necessary to know 
the over-all multiplicity m(E, x). By making an assumption that m(E, x) is propor- 
tional to E“, where a is constant, it is possible to obtain m(E, x) of the ionizing and 


=. 
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neutron components by using the observations made by H.V. Neher and J.A. Simpson 
et al., respectively (see App. II). The results obtained are shown in Figs. 1 and 2, and 
in Table I. Two different m(E, 1030) of the ionizing component are tabulated in 
Table I. These over-all multiplicities are those of the two extreme cases estimated 
on the basis of assumption (see App. II). The full lines in Figs. 1 and 2 denote the 
over-all multiplicities obtained directly from the observations without any assumption 
and E, in Table I is the critical energy at which the distribution of m(E, x) changes. 


E (Bev) 
Fig. 2 Over-all multiplicity of the neut- 
1 10 10° ron component. 


E (Bev) 


Fig. 1 Over-all multiplicity of the ioni’ 
zing component. 
Table I m/(E, x) 


x (gcm-*) m(E, x) (ECE,) E,(Bev.) | m(E, x) (ExE;) 


= 0 1 1 
I E \0. 
a 150 See Fig. 1. 10 6.1-( 
a | 2 
a | é ye3 
S 300 7.4+10-2-E18 8.2 3.0 ( = 
oo ; : 6 
4 I 1.3-10-4-E38 15 2.8+10-1- ex 
2 | 1030  Lceelteeer as 2S S Fok 
sa I 1.3-107-E*5 23.1 | 9 5-10-1.( e 
e.4 305 See Fig, 2. 4 50+ £0.28 
Bae Mie ama 
36 680 * See Fig. 2. 7 8.6- E08 


- ; 
* (This over-all multiplicity was oftained from the latitude curve of intensity nor- 
malized to unity at the equator.) 


In Figs. 3 and 4, the functions of n(E, x)=m(E, x)t(E) of the ionizing and 
neutron components are shown respectively and normalized to unity at E=10 Bev. to 


facilitate their mutual comparison, 
e 
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Fig. 3 Relative intensity spectrum, 
n(E, x)=m(E, x)-i(E), of 
ionizing component. 
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Fig. 5 Latitude dependence of y of the 
ionizing component by the varia- 


tion of the_staticzelectric, field,j 
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Fig. 4 Relative intensity spectrum, 
n(E, x)=m(E, x)-i(E), of the 
neutron component. 
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Fig. 6 Latitude depeddence of y of the 


neutron component by the varia- 
tion of static electric field. 
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By using these over-all multiplicities and the primary spectrum given by Eq. (6), 
it is possible to calculate y(A, x, 4E, 0). Figs. 5 and 6 show the latitude curves of y of 
the ionizing and neutron components at various altitudes in case 4E=0.1 Bev. In 
these figures, the subscripts 7 and n of y represent the ionizing and neutron ae 
nents respectively. 

The altitude curves of y at various latitudes are obtained by drawing the 
smoothed curves through the calculated points and shown in Figs. 7 and 8 In 
Fig. 7, the curves of the ionizing component are normalized to unity at x=0. 


yn(Q, %, 0.1, 0)/y¢ (2, 0, 0.1, 0) 


v0, x, 0.1, 0)/y« Q, 0, O,4:, 0) 


10°" 


5-10? 10° 
Atmospheric depth x (g cm-) 
Fig. 8 Altitude dependence or y of the 
neutron component by the varia- 
tion of the static electric field. 


0 5-102 10° 
Atmospheric depth x (g cm-?) 


Fig. 7 Altitude dependence of y of the 
ionizing component by the varia- 
tion of the static electric field. 


The variation of the neutron component is produced by the variation of the 
primary ionizing particles which is caused by the variation of the electric field, 
so that the curves of the neutron component are normalized so as to satisfy the 
condition that the variation of the primary particles becomes just to unity. J.A. 
Simpson et al.“*) reported that, in the atmospheric depth greater than 600 gcm~’, 
there were no differences in the latitude dependence of the neutron intensity. If so, 
the variation y of the neutron component must be independent on altitude in such a 
region. On the contrary, Fig. 8 obtained by drawing the smoothed curves through 
the calculated points shows that y decreases slowly with the increase of x in such a 
region. This discrepancy between the theoretical and experimentally expected curves 
might be due to a lack of the calculated points, 
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g 3. The Cosmic Ray Intensity Variation 
by the Variation of Cut-Off Energy 


As discussed in § 1, the variation of cut-off energy is another cause which 
produces the cosmic ray intensity variation, and cannot be uniquely determined unless 
the state of the static electric and magnetic field surrounding the earth is determined. 
Even if the state is determined, it is still difficult to obtain the variation of cut-off 
energy. In App. I, this problem is treated in a simple case when the static electric 
and magnetic field is symmetrical to the earth’s magnetic dipole axis. 

The latitude curve of y cannot be obtained when the variation of cut-off energy 
is not known as a function of latitude, but it is possible to obtain the altitude curve 
of y if only N(A, x) can be known from the observations. If the cut-off energy FE, 
changes to F,+6E, by the influence of the static electric and magnetic field, y in Eq. 
(5) becomes as follows, 

E)+6E) 
| m(E, x)i(E)dE 
¥(A, x, 0, 0F,) = wee oe po ae (9) 

Figs. 9 and 10 represent the variation y at various altitudes and latitudes in case 
OF, =0.1 Bev. Care must be taken that, as mentioned above, these curves do not 
represent those of the latitude effect of y, different from those in Figs. 5 and 6. At 
high latitudes, y=0 as shown in Figs. 9 and 10. This phenomenon is partly due to the 
fact that there are no primary cosmic rays in the lower energy region other than the 
general cut-off energy E, and partly to the fact that at the low atmosphere there is 
the knee of the latitude effect of cosmic ray intensity produced by atmospheric 
absorption. But if 0d, becomes very large under the influence of some mechanism 
even in such high latitudes, the variation y occurs as well as at low latitudes. 


2 


yiQr, ¥, 0, 0.1) (%) 


nO, x, 0, 0.1) (%) 
— 


0 30 60 90 0 30 60 90 
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Fig. 9 a y, of the ionizing component Fig. 10 Variation, y, of the neutron component 
y the variation of the cut-off energy. by the variation of the cut-off energy. 
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Figs. 11 and 12 show the altitude curves of y produced by the variation of cut- 
off energy. The normalization is quite similar to that adopted in case of Figs. 7 and 
8. As seen from Figs. 7, 8, 11 and 12, the altitude curves, in case of the variation of 
cut-off energy, decrease more rapidly with the increase of x, compared with those in’ 
case of the variation of electric field. Such a difference is caused by the fact that, in 
the former case, only the low energy particles are excluded; in the latter case, on 
the contrary, all the particles are influenced by the electric field. 
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Fig. 12 Altitude dependence of y of the neutron 
component by the variation of the cut-off 
energy. 
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Fig. 11 Altitude dependence of y of the ionizing 
component by the variation of the cut-off 


energy. 
§ 4, Comparison with Observations and Discussion 


Based on the calculations in § 2 and $ 3, we consider which variation of either 
the electric field or the cut-off energy is best suited to explain the world-wide varia- 
tion of cosmic ray intensity (not dependent on local time) which has three character- 
istics as cited in the Introduction. y 


[1] The variation is world-wide. 
In case of the electric field:--This can be explained if only the potential on the 


earth varies. Moreover, as pointed out by H.V. Neher,“ the fluctuation in intensity 
must be the same all over the knee of the latitude curve of intensity. This pheno- 
menon can be also explained, independent of the formation of the knee, as shown in 


Figs. 5 and 6. 
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In case of the cut-off energy :—At low latitudes, the changes occur without any 
condition if the state of the static electric and magnetic field varies, but some difficulty 
arises in high latitudes. As pointed out by various authors,“ ©“ the change cannot 
be produced over the knee by the variation of the cut-off energy. They are right if 
the static electric and magnetic field is symmetrical to the earth’s magnetic dipole 
axis, since such a field, however large, does not produce any variation of cut-off 
energy near the magnetic poles and then does not satisfy the condition that the 
variation must be the same all over the knee. On the contrary, if we consider 
the field not to be symmetrical to the earth’s magnetic dipole axis, such a field 
has the possiblity of satisfying the above mentioned condition. The reason is that 
this field, a simple example of which is the solar magnetic field, has an agency to cut 
off the low energy region of the primary spectrum,@” so that, if such a field is 
produced or changed by some mechanism, the general cut-off energy E, is changed 
and then the intensity variation is produced and is the same all over the knee. But 
a difficulty arises even in such a case; at the normal state, the knee at low altitude 
is caused by the atmospheric absorption and not by the general cut-off energy E,, so 
that, in expectation of the intensity decrease over the knee at such a low altitude, 
the general cut-off energy must become greater than the cut-off energy FE, at the 
latitude where the knee at low altitude begins. If such a variation of E, is produced, 
the latitudes, at which the knees at various altitudes begin, must become constant. 
Such a phenomenon has not yet been reported. Moreover, we cannot expected any 
intensity increase from the normal state over the knee at low altitude by such a field, 
since the decrease of E, gives no influence upon the intensity over the knee at low 
altitude. 

As discussed above, it seems to be difficult to explain the phenomenon [I] by 
the variation of cut-off energy, although there is some slight possibility of doing so. 

[II] The latitude dependence of the variation 

Up to date, not much has yet been done about the latitude effect of the world- 
wide variation at sea level. This is due to the following difficulties: (1) The observa- 
tion stations are not located at the same altitude, (2) each station is located at differ- 
ent longitude and (3) the type of the cosmic ray meter is different at each station. 
As we know nothing about these corrections for the variation, the true latitude 
dependence of y at sea level cannot be obtained. In these corrections, (1) gives the 
most important influence upon the latitude dependence of y. 

In case of the electric field:—Keeping the above in mind, we discuss the 
latitude dependence of y, based on the observed results obtained by several authors.@» 
@O.@2) The relative values of y, at each station, normalized to that at Cheltenham, are 
shown in Table II. As seen from this table, at some stations, there are some differ- 
ences between the values obtained by these authors. For example, the value at 
Teoloyucan obtained by S.E. Forbush is exceedingly larger than that obtained by 
S. Yoshida and Y. Kamiya, and another example is that the value at Freiburg is 
extremely small compared to the values at other stations which are located at: the 
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same altitude and nearly the same latitude as Freiburg. On the other hand, it is a 
remarkable point that at Huancayo there is quite a good agreement between them. 
To compare these data with the theoretical latitude curve of yx: in Fig. 5, it is neces- 
sary to make the altitude correction to the values at several stations (i.e. Huancayo, 
Teoloyucan, Canberra and Hafelekar). 


Table II Relative value of y; at each station normalized to that at Cheltenham and not 
correted the altitude effect. 


Geomagn.| Height above Relative value of ys obtained by 


Latitude | sea level (m) - ee IM. Wada@2* 


No. Station 


S.E. W.H. 
Forbush | Fonger® § 


1 Godhayn 80°N 1.124.13 0.97 
2 Bergen 61 N war ley i aateepe rie, 
3 Manchester 57 N eee Soeeeee Rak rmeN cre a 
/ 4 Boston 54 N 0.97 
5 Cheltenham 50 N 1.00 Lo | Yio.) ole 
6 Christchurch 49 § | eeern 0.95 
7 Friedrichshafen 49 N | 
8 Freiburg 49 N 0.77 
9 Hafelekar 48 N 1.334.083 ele Baie UE 
10 Canberra 45 S 0.944+.15 
11 Teoloyucan 30 N 1.18-+.08 1.42 
12 Tokyo 25 N 0.85.05 0.93 
13 East China Sea| 22 N 0.87 
14 Huancayo ne 0.894.07 0.92 0.90 0.91 


* The values in this column were obtained from the data of cosmic ray storms. 
§ The values in this column were obtained from the data of 27-day variations. 


Let ¢;(A, x/—x, 4E,0) denote the altitude correction factor to alter y; at x to 
that at x’, and the following relation holds 


(A, 1030, JE, 0)=c;(A, 1030—x, 4E, 0) y;(4, x, JE, 0). (j=7 or n) ----(10) 


c; is obtained theoretically by using the altitude curve of y; as shown in Figs. 7 and 
8, and is listed in the 4th and 5th columns of Table III for each station in case j=7 
and 4E=0.1 Bev. 

Using the altitude correction factor c,’, we can obtain the experimental relative 
value of y; at sea level. The results are shown in Fig. 13 together with the latitude curve 
of y,7 (4, 1030, 0.1, 0) obtained theoretically. As seen from this figure, there is a fairly 
good agreement between the observations and the theoretical curve, although the 
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Table IIL Altitude correction factor ¢; (>, 1030—-x, 0.1, 0) 


@1) 


Station aaa ele A ot : ct ” Observed value of ¢; 
Boulder 49°N 1600 0.77 0.83 

Climax 48 N 3500 0.60 0.70 

Hafelekar 48 N 2300 0.69 0.77 | 0.71 (Cheltenham/Hafelekar) 
Bagnéres 46 N 550 0.94 0.95 

Canberra 45S 800 0.87 0.92 | 1.06 (Cheltenham/Canberra) 
Teoloyucan 30 N 2285 0.77 0.87 | 0.725 (Tokyo/Teoloyucan) 
Huancayo : 1S 3350 0.80 0.96 


* J and I denote the correction factors obtained by using the value y:7 (2, 1030, 0.1, 0) 
and y:// (2, 1030, 0.1, 0) respectively, as shown in Figs. 5 and 7. 
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Fig. 13 Latitude dependence of »; at sea level, nor- 
malized to that at Cheltenham. Numbers 
denote Stations in Table III. 
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differences between them are 
large at some stations (for 
example, Freiburg and Canberra). 
The experimental latitude de- 
pendence of y; at sea level, 
corrected by c,“, also has a fair- 
ly good agreement with the 
theoretical one, but is not shown 
here. As pointed out previously, 
there is quite a good agreement 
between the experimental values 
of y, at Huancayo, obtained by 
various authors. The value 
without the altitude correction 
is equal to 0.89-40.07 which was 
obtained by S. Yoshida and Y. 
Kamiya. On the other hand, the 
theoretical value is obtained 
from Figs. 5 and 7 and equal to 
0.82 in case y,’ is adopted and 


0.78 in case y,;“ is adopted. The theoretical value in the former case is better than 
that in the latter, and is in accordance with the observed value within the experi- 
mental error. From Figs. 5 and 6, we see that the latitude dependence of y at high 


altitude satisfies qualitatively the condition that the changes must increase with the 
increase of latitude. 


S. Yoshida and Y. Kamiya@” obtained the altitude correction factors at several 
stations from the observation by comparing the values of y; at two stations which 
were located at different altitudes respectively and nearly the same latitude each 
other. The results are tabulated in the last column of Table III. We see that there 
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are fairly good agreements between the theoretical and experimental altitude correc- 
tion factors except that at Canberra. Although it is difficult to interpret that the 
experimental value at Canberra becomes larger than unity, this might be due to the 
scarcity of observations or the difference of cosmic ray meter compared with that at 
other stations. 

As seen from the above comparisons of the theoretical values with the obser- 
vations, the latitude dependence of the world-wide variation is explained by the 
variation of the static electric field. 

In case of the cut-off energy :—As pointed out in § 3, the latitude dependence 
of the variation cannot be obtained unless the state of the electro-magnetic field is 
determined, so that the comparison of the theoretical curve with the observation 


cannot be done. 

{III} The altitude dependence of the variation 

Recently, the world-wide variations of the ionizing and neutron components at 
various altitudes have been observed by various authors.©>@3.@) H.V. Neher pointed 
out that the ratio of the variation at high altitude to the variation at sea level is 
always constant. Such ratios obtained experimentally are reproduced in the 3rd 
column of Table IV. The theoretical ratios expected from the variation of the static 
electric field are in the 4th column where I and II represent the values obtained by 
using y,/(A, 1030, 0.1,0) and y,/{A, 1030, 0.1, 0), respectively. In case of the variation 
of cut-off energy, it is possible to obtain the ratio of the variation at the same 


Table IV Altitude dependence of y 


| 
Censeeed Value expected from| Value expected 


yi, %) 2 920, *) rcaicar electric field from cut-off Remarks 
yi(56, 50) : (56, 1030) veh Ftd. a6: 1 (3) 
Bismarck 
yi(45, 50): y(1, 700 ) Cok 4.6:1 aati: 1 50:1 (3) 
Chicago Huancayo 
; ; ; (3) 
yi(43, 50) : yi(41, 880 ) 4<1 Ze L 7A) Fab Maal ike cats arise 
Be Waeeie* <Sies Wiles obtained at the 
vee time of the 
magnetic storm 
¥yn(48, 680) : y:(50, 1030) | 4.0:1 See e oe key © 00% 1 (6) 
Climax Cheltenham 
yn(48, 680) : yi(49, 1000) | 5.5:1 (6) 
Climax Freiburg 
Yyn(48, 680) : yi(1, 700) Ais A 4,821 A, Oval 50.: 1 (6) 
Climax Huancayo 


- Numbers in the remarks column denote those of References quoted. 


latitude, but not possible to obtain the ratio if stations are located at different lati- 
tudes, because, as stated in § 3, the latitude curve of y» cannot be obtained unless 
the state of the electro-magnetic field is determined. But, if it is assumed that the 
intensity variation caused by the variation of cut-off energy shows no latitude effect 
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at ground level as the world-wide variation actually observed at the ground shows no 
remarkable latitude effect, it is possible to obtain the ratio between two stations which 
are located at different latitudes. The results are shown in the 5th column of Table 
IV. As stated above in [I], the variation of cut-off energy could possibly produce the 
intensity variation at high latitude, but here no such ratio is obtained. 

As seen from this Table, the observed ratios show better agreements with the 
theoretical ones expected from the variation of electric field than that expected from 
the variation of cut-off energy. Moreover, in case of the electric field, the values in 
case I are better than those in case JJ. But, in point of detail, there are some 
discrepancies observable between the theoretical and experimental values as shown in 
the 2nd and 3rd rows of Table IV. These discrepancies seem to increase with the 
decrease of latitude. ‘ 

In the 4th, 5th and 6th rows, the ratios between the neutron variation at 
Climax and the ionizing variations at several stations are shown. The observed ratio 
adopted by W.H. Fonger in order to compare with his theoretical calculation is that 
cited in the 5th row. But, as pointed out previously and seen from Table II and 
Fig. 13, the variation of the ionizing component at Freiburg is very small, compared: 
with the variations at other stations which are located at the same latitude as 
Freiburg, so that the ratio in the 5th row is much larger than that obtained by using 
the variation at the other stations. Then, in order to compare with the theoretical 
value, it is more reasonable to adopt the values given in the 4th and 6th rows, than 
to adopt the value given in the 5th row. As for the ratio of the neutron component 
at Climax to the ionizing component at Cheltenham or Huancayo, there is quite a good 
agreement between the experimental and theoretical values in case I. 

As stated in § 2, the variation of the neutron component must be independent 
on altitude in the atmospheric depth greater than 600 gcm~, if the observation of 
the neutron intensity made by J.A. Simpson et al.“*) is correct. Therefore, to confirm 
this fact, it is most desirable to make the comparisons between the variation of the 
neutron component at Climax and those at the other stations. 

As seen from the above various considerations and discussions, it is possible to 
explain all the phenomena of the world-wide variation by the variation of the earth’s 
potential, although in point of detail there are some discrepancies observable between 
the experimental and theoretical values. On the contrary, the variation of cut-off 
energy cannot explain the altitude dependence even though it may possibly satisfy 
the phenomenon [I]. 

[IV] The large variation of cosmic ray intensity 

The large decrease of 6% or more in cosmic ray intensity at sea level are often 
observable at the time of the cosmic ray storm.@ In such cases, the altitude and 
latitude dependences of the variation calculated in § 2 do not hold, because, as pointed 
out in § 1, they depend on the magnitude of the variation of earth’s potential.“ We 
must therefore calculate them when the variation of earth’s potential is large. Figs. 14 
and 15 represent the latitude dependences of y and Figs. 16 and 17, the altitude 
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dependences of y, in case 4E=1 Bev. The method of the normalization of the curves 
in Figs. 16 and 17 is the same as that adopted in Figs. 7 and 8. From these figures, 
we see that the altitude and latitude dependences become smaller if 4E becomes 
larger. Moreover, to make this feature more clear, Figs. 18 and 19 are shown, 
representing the variations of the ionizing and neutron components at various altitudes 
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Fig. 15 Latitude dependence of y of the 
F neutron component. 
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Fig. 14 Latitude dependence of y of the 
ionizing component. 
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as a function of JE, in case 458°. As seen from these figures, the dependence of 
y on 4E becomes smaller with the increase of x, and it can be concluded that y is 
approximately proportional to 4E at low altitudes and latitudes. But, at high altitudes 
and latitudes, it depends exceedingly on the magnitude of 4E. Then it is necessary 


— 
(om) 


10 


ou 


yn(i>58°, x, AE, 0)/yn (458°, ¥, 0. 1, 0) 


yiQ>58°, x, AE, 0)/94 (.>58°, x, 0-1, 0) 
on 
= 
en & 


0 0 
0 2 1 0 0.5 1 
AE (Bev) : AE (Bev) 
Fig. 18 AE dependence of y of the ionizing Fig. 19 AE dependence of y of the neutron 


component. component. 


to consider such an influence upon the variation in order to compare the experimental 
altitude dependence of y with the theoretical one. 


§ 5. The Diurnal Variation of Cosmic Ray Intensity 


The diurnal variation of cosmic ray intensity is also a world-wide phenomenon 
at sea level, whose difference from the world-wide variation discussed in the previous 
sections is that the former depends on local time. Here, an examination is to be 
made of the possibility that the diurnal variation can be explained by an electric 
field@® although the origin of such an electric field cannot be pointed out at present; 
and a qualitative comparison of the diurnal variation expected from such an electric 
field hypothesis is to be made with that expected from the static magnetic field 
hypothesis.‘” 

The electric field derived from a potential can no longer explain the variation 
which depends on local time, because the variation is produced simultaneously all over 
the world as long as the electric field is derived from a potential.¢? Then, to explain 
such a variation, it is necessary to consider an electric field not ‘derived from a 
potential. In such a general electro-magnetic field as well as the static electric and 
magnetic field, the intensity of cosmic rays, whose momentum are P at any point, is 
also expressed by , 


$oP*7). 08, Oya ee tala as (11) 


where D=SN, with dv=dxdybzdP,6PyBP, , +-+0-2 +s (12) 
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and dN is the number of particles in the volume dv. Even in such a general field, 
the extended Liouville’s theorem®>©® holds and then D in Eq. (11) is constant as well 
as in the static electric and magnetic field.“ But, the difference in both cases is 
that, in the latter case, P is determined if only the electric potential at any point is 
known; in the former case, on the contrary, P cannot be determined unless the 
trajectory of charged particle can be obtained in the general electro-magnetic field. 
Keeping this point in mind and assuming that the primary particles consist of protons 
only, we are able to obtain also the vertical cosmic ray intensity at geomagnetic 
latitude 4° and atmospheric depth xgcm-* in the general field, using the method 
similar to that adopted in the static field. 


NU, x, AEE, 2,1), 8E}=\ | m(E, 2)[1-L{B+ AEE, 4, 1), AEE, 2,9) 
a a 
xt{E+4E(E,1,d}dE, ----(13) 


where N,6E,, m,L and i are the same notations as those in Eq. (3), and 4E(E,A,2) 
is the energy loss which the primary particle, having its kinetic energy E at 4° on 
the earth, undergoes during its passage through the electro-magnetic field, and is a 
function of E, 4 and local time ¢, since 4E depends upon the particle’s trajectory 
which can be determined by parameters EF, 4 and ¢. Although it is meaningless to 
try to determine whether such a variation is caused by the electric field or the 
magnetic field, we designate such a variation for the variation caused by the general 
electric field in distinction from the variation caused by the static electric field or by 
the static magnetic field, since the agency of the electro-magnetic field upon the intensity 
variation is similar to that of the static electric field as seen from Eqs. (8) and (13). 
Fig. 20 isa schematic representation of 
deformations of the primary spectrum which 
will be caused by the general electric field. 
The distinction between the intensity variation 
caused by the static electric field and that 
caused by the general electric field is as 
follows: In the former case, the relative 
contribution of the primary particles, each 


having a different energy, to the intensity 


Fig. 20 Variation of primary - — variation is uniquely determined (see Eq. (8) 
spectrum produced by the electric A 
field. JG -! and Curve B in Fig. 20). In the latter case, 
1. Curve A represents the spect- on the contrary, the relative contribution 


rum at the normal state. : ‘ 
2. Curve B represents the varia- varies with local time (see Eq. (13) and 


tion of the spectrum from the Fe - 
atl ieltas wradeced hip. the Curves C and C’ in Fig. 20). oso me 
static electric field. altitude dependence of diurnal variation is 

3. Curve C and C’ represent the 
variations of the spectrum. pro- 
duced by the general electric tion expected from the static electric field 
field at the times t and t’ res- 


pectively. 


different from that of the world-wide varia- 


and cannot be obtained theoretically unless 
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the state of the electo-magnetic field is determined. 

Comparing the diurnal variation of neutron component at Climax with that of 
ionizing component at Freiburg, W.H. Fonger“ showed that the ratio between the 
amplitudes of variation at the two stations is nearly equal to the ratio of the world- 
wide variation (not dependent on local time) expected from the static electric field 
hypothesis (here-after referred to as Phenomenon b). If we substitute the altitude 
dependence curve obtained from the static electric field for the true altitude depend- 
ence curve of diurnal variation as a rough approximation because of the existence of 
Phenomenon b, the amplitude of diurnal variation of the ionizing component at high 
latitude will be expected to be about 10 times that of sea level at the top of the 
atmosphere and about 7 times at *=50 gcm~’. 

The static electric and magnetic field also has the possibiblity of producing the 
diurnal variation, if the field is not symmetrical to the earth’s magnetic dipole axis. 
There is the solar magnetic field theory adovocated by M.S. Vallarta,“” as the repre- 
sentative example. The altitude dependence of diurnal variation expected from this 
theory is smaller than that expected from the variation of cut-off energy, because, as 
pointed out in § 1, the diurnal variation is caused by the variation of primary spect- 
rum at intermediate energy region and not by the variation of cut-off energy. The 
distinction between the general electric field hypothesis and the solar magnetic field 
theory is as follows: The altitude dependence of diurnal variation expected from 
the former is smaller than that expected from the latter,“ since, in the latter case, 
only the primary particles in the low energy region are influenced; in the former 
case, on the contrary, all the particles in all of the energy region are influenced. 

As above, the diurnal variation could possibly be explained by the electric field 
not derived from potential, even though it cannot indicate the mechanism producing 
such an electric field. Recently, the diurnal variation at high altitude was observed 
by various authors“ in order to confirm the existence of solar magnetic field. But, 
none of them has reported the ratio between the amplitude of diurnal variation at 
high altitude and that at sea level. It is most important to examine the altitude 
dependence from the experimental and theoretical points of view in order to determine 
the cause of the diurnal variation. 


Concluding Remarks 


It is concluded that the world-wide variation of cosmic ray intensity not depend- 
ent on local time can be explained by the variation in the distribution of the primary 
cosmic ray intensity spectrum which is caused by the variation of the static electric 
field, but not explained by the variation of cut-off energy which is caused by the 
variation of the electro-magnetic field. 

The diurnal variation of cosmic ray intensity could possibly be explained by an 
electric field not derived from potential. The existence of Phenomenon b seems to be 
favourable to the general electric field hypothesis. 
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Appendix I The Variation of Cut-Off Energy 
by the Static Electric and Magnetic Field 


The relativisitic equation of the charged particle’s motion in the electro-magnetic 
field is 


d " e 
a (me) =e + (oe (A, 1) 


General solution of Eq. (A,1) cannot be obtained analytically even in case of the 
motion in the earth’s dipole magnetic field. Then it is extremely difficult to solve 
Eq. (A, 1) in case the static electric and magnetic field surrounding the earth exists in 
addition to the earth’s magnetic field. But, in one special case, it is possible to obtain 
the Stérmer cone“ and then to obtain the varia- 
tion of cut-off energy caused by the static electric 
and magnetic field. As the coordinate system, here 
are the spherical coordinates (R, 4, v) adopted, as 
shown in Fig. 21, and the z-axis coinciding with the 
earth’s magnetic dipole axis and directed towards 
the North Pole. 

If the static electric and magnetic field sur- 
rounding the earth is symmetrical to z-axis, the 
following conditions hold; 


E,=0 and H,=0. -:----- (A, 2) 


Fig. 21 Coordinate system 


Then, in such a case, the ¢g-component of Eq. (A, 1) 
becomes 


d — é eoeeovr ewe eeoe 
a (m D)» a (vuH],.- (A, 3) 
In Eqs. (A, 1) and (A,3), the mass m is not constant because of the existence of the 


electric field. But, quite similar to the case that only the static magnetic field 
exists, Eq. (A, 3) can be integrated, and the solution is 


m(Reos A) Se = ta — (Roos) Ay +7"| Peat Cire Pr (A, 4) 
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where J” is an integral constant, and A, is the y-component of vector potential and 
in this case is given by 

A, = —Moos*_ as, cee (A, 5) 
The Ist term on the right hand side of Eq. (A,5) is due to the earth’s dipole 
magnetic field, M denoting the magnetic dipole moment, and —A% to the outer field. 
Let 0 denote the angle between the particle’s orbit and the meridian plane, and the 
next relation is obtained by using Eq: (A, 4). 


. » Rcosi dy __e {Mcosh Oy ee 
Sint ven imal Ie hls s(RA+ Regsi | eP 


Then, the motion must satisfy the following condition : 


e {Mcosh Pl nee 
moc | + A(R A) Re eost f S* (A2 


Next, we consider the extended Stérmer transformation®“® as follows, 


|e| M ise pig. vsaheenMeres 
ane and Re=lr, etc., (A, 8) 
and 
(mvc) =(mve)[1—L{ Ep, €36(R, AYP, vere cree eee (A, 9) 


where ZL is Liouville’s effect and given by Eq. (4), E is the kinetic energy correspond- 
ing to the value of mvc, and subscript , denotes the value on the earth. 06¢(R, 4) is 
the difference in electric potential relative to the earth’s potential at any point (R, A). 
Appling these transformations to Eqs. (A, 6) and (A, 7), the following relations can be 
obtained ; 


1 is Be : 

sind=q— pn “reel eka) ep ee 
cos A 271° 

Fresaue rd poy t pp Avcer, ah j |= sreseses (A7Y 


where 27, is the integral constant. From Eq. (A,7/, we can obtain the following 
relation ; 


A 
reosi( “S +1l+a— ~8)=2r=reosi( 5414048), tenes (A, 10) 
where 
alr, &y ro) 7 AX(tr, a), 
: (A, 11) 
B (7, A, Yo) ad 1—[1—L{E,, ed (ir, ay" e 


a and # are the terms due to the outer magnetic and electric fields respectively, and 
are functions of 7, A, and 7% (which is the particle’s energy in St6érmer unit@® cor- 
responding to FE) or (mvc))). Eq. (A, 10) depends upon the particle’s energy, different 
from the case where only the earth’s magnetic field exists. Hence, the Stérmer 
transformation loses its usefulness in this case. 
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As the influence of the outer magnetic field upon the cut-off energy has been 
examined by various authors,“ here that of the electric field only is considered. In 
this case, Eq. (A, 10) becomes as follows, 


reosi{ S41 8)=2r=roosa( S414 8). tees eees (A, 12) 


Fig. 22 is a schematic representation of Eq. (A,12) in case 4=0. The full lines 


ui 


= represent the curves n=4(4+ 1) 

and n=$(4- 1), which equations 
} aphiddes denote the non-existence of the electric 
region field and are tenable independent of 
the particle’s energy under considera- 
tion. When the electric field exists, 
these curves cease to be independent 
of the particle’s energy. The dotted 
lines in Fig. 22 represent an example 
of the deviations from the normal 
state produced under the influence of 
the electric field, in case of the 


particle’s energy becoming 7” Stérmer 
forbidden 


region (7)<1) on the earth. To examine the 


@=0 variation of cut-off energy in such a 
Fig. 22 yi—7 diagram case, it is necessary to obtain the 
minimum value of the following curve; 


2n=rt1—Ble, 0, rd}. pn eine owe (A, 13) 


Generally, this problem cannot be solved unless the form of 8 is determined. But an 
analysis can proceed when the disturbance by the electric is small. Under such 
a condition, 8 in Eq. (A,11) is approximately given by 


B(r, 0, 10)~ Bap (rs 0, 10) = L{Eo, e86(r, 0, 79)} 


oA) er 
(Bot me? — (me)? sa : ; 


where E, is the particle’s energy in Bev., corresponding to 7 Stormer. The above 
equation holds if only E,>ed¢(7, 0, 7%). As the minimum value of 7; in Eq. (A, 13) is 1 
at r=1 in case B=0, it can be assumed that the minimum value of 7; will be given 


by 


(71) mim =1+471(%, 8), (471<1) 
at 


Y =1+47(%, B). (4r<1) 


Let @ be further developed in the power series at r=1: 
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/ A VE 0, i) 2 $ . 
A(r, 0,7») 81, 0,79) + 2.0) (r—1) + 2, Oto) yay fo, ++ (A, 16) 


where ((1, 0, 7) is the value at y=1 and §’(1, 0, 7) is that of the 1st derivative of B 
with respect to 7. Using Eqs. (A, 13), (A, 15) and (A, 16), and neglecting the higher 
order of dv as small, we can obtain the following relations : 


1 {B (1, 0, 7o)+6'0, 0, Yo}? 
Ari(o Byes ->|¢ CL 0, Yo) +3970, 7o.)—b"(0,0, I. 


= BQ, 0, 7%))+6'(0, 0, To) 
Avr, 8) F= 28, 0, 74) BUs 0, 74) 


ies (AD 


Next, we examine the variation of cut-off energy produced by such a variation of the 
minimum value of 7;. In case only the electric field exists, Eq. (A, 6)’ becomes as 


follows, 


, 


in@=—1—.( 2084 __ 271 ) 7 ae 

sin#= (i—D)' ae 4 (A, 18) 
Introducing, into the above, the values 0=0, r=7%, 11=(71)min» and Z=0 (as the 
potential difference 6¢ is zero on the earth), the following relation is obtained, 


cosk=plreisind+V Ge Sin OF FBI Tm +e+++++ (A,19) 


A,, satisfying Eq. (A, 19), denotes the latitude at which 7) Stérmer, being the particle’s 
energy on the earth, the particle entering into the earth from the given direction 8@, 
exactly becomes the cut-off energy. Incase the electric field does not exist, let 4) be 
the latitude at which 7) St6érmer becomes the cut-off energy, and dy is given by 


cosky=4 [r'sin 6+ (r,2sin 02 +875]. tereeeeeeees (A, 20) 


From Eqs. (A, 19) and (A, 20), we get 


2% ° Ari “ees Ow vedic wan (A, 21) 


here a 
iif vf V (ro? sin? +875’ 


cos 4, ~cosdy+7, ! 


in which 4;, is given by Eq. (A,17). Moreover, in case the electric field does not 
exist, let 7, Stérmer be the cut-off energy at the latitude 2,, we get the following 
relation,@” 


71=COSAy[1+. UM cos*asinO]-1, sven nn enedee (A, 22) - 


Introducing Eq. (A, 21) into the above, the following relation holds between 7, and 7%, 


coskg 2 rocosdysind+1 


71 ~7 (1+), 
where 3 ( 2 3 7sin@ Jp. 4's ce (A, 23) 


If the energies 7; and 7) Stérmers in Eq. (A, 23) are converted into the energies E, 
and E, in Bev. respectively, Eq. (A, 23) becomes 
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Bong Bete te 
Ww ere - (Eo mye? —(mne?)? 
¢= Eotmec -(20),; Cres eees (A, 24) 
and: - | 
ca (nr) [(Eo+ moc??? — (mac*)*1*, 


in which R, is the earth’s radius. From Eqs. (A,21) and (A, 24), we see that the 
particle’s energy E,, which is the cut-off energy at 4=A, and 0=0 in case of the non- 
existence of the electric field, becomes the cut-off energy at A=A, and 0=0 in case of 
the existence of the electric field, and that the relation between E, and FE, (which is 
the cut-off energy at A=A4, and @=0 in case of the non-existence of the electric field) 
is given by Eq. (A, 24). 

To make the physical meaning of these equations more clear, the simplification 
of these equations is carried out. In Eq. (A,17), 47; is a function of 8, 6’ and ~”, 
and then cannot be determined unless the space distribution of the electric potential 
is known. But, if we assume that the potential does not vary rapidly near the point 
where r=1 and 4=0 so as to satisfy the following relations, 


B’aR<1 and b’<P<l, 
and then, neglecting the higher order of f, #’ and ~”, Eq. (A, 17) becomes as follows, 


“1: (Eo +myc*) 
2 (£ot+moc*)?—(ame’)? 


4n~—5 BU, 0, %)=— -e0¢(1, 0, 10). ie) -(A, 25) 


Using the above equation, in case of the vertical incidence (@=0), Eq. (A, 24) becomes 


as follows, 
F,=E,—e6¢(1, 0,7). Bev.  — sveeveeeveee (A, 26) 


As seen from the above, the variation of the cut-off energy produced by the static 
electric field is approximately determined only by the earth’s potential, and is 
—ed¢(1, 0, 7) in case of the vertical incidence. When the potential at the point where 
y=1 and A4=0 is higher than that of the earth (d¢>0), the cut-off energy becomes 
larger than that of the normal state and then, as discussed in § 3, the cosmic ray 
intensity is decreased by such a variation. But, as seen from Fig. 9, such a decrease 
is exceedingly small even in the high atmosphere if ed¢<10-? Bev. As pointed out 
by H.V. Neher, there exists the general cut-off energy E,(=0.8 Bev. for proton) 
which corresponds to about 0.16 Stérmer. The distance from the earth’s center to 
the point where v=1 becomes nearly 6 times the earth’s radius even if we consider 
the lowest energy 0.16 Stérmer. Therefore, if the following relation is satisfied over 
the region less than 6 times the earth’s radius, 


e0¢2107 Bev., 


the cosmic ray intensity variation is not produced over all the latitude by the 


variation of cut-off energy. 
In conclusion, the variation of cut-off energy can be solved analytically in case 


the static electric and magnetic field is symmetrical to the earth’s magnetic dipole 
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axis. The variation of cut-off energy caused by the static electric field can be 
determined approximately by the potential difference at the point where y=1 and 
A=0, relative to the earth’s potential, as long as the disturbance by the electric field 


is small. 
Appendix II Determination of Over-all Multiplicity 
(1) In case of the ionizing component 
From Eq. (7), we can obtain the following relation 


___1_.dNQ,*), di 
Oe EY dh. are 


H.V. Neher“ obtained m(E, x) at low energy region, using the primary spectrum 
in Eq. (6) and the vertical intensity of the ionizing component at various altitudes. 
The results are reproduced in Fig. 1 as the full lines. In the high energy region, 
although m(E, x) cannot be obtained directly from observation, it is possible to obtain 
it approximately under the following considerations: As shown in Fig. 1, in case x 
is small, the full line has a critical energy E,(x), in the energy region larger than 
which m(E, x) changes more slowly compared with that in the low energy region. 
In such a high energy region where E>E,(x), we assume that m(E, x) is expressed 
by the following formula, — 


m(E,x)=m42)| aos | : sete ee eeeeee (A, 28) 
E{x) 
where (x) is the value of m(E, x) in case E=E,(x). Introducing Eq. (A, 28) into 
Eq. (7), we obtain the following relation, 

NQax)=5~ mA) acy | tats peel ial (A, 29) 

Ec) Ex) 

where 4, is the latitude at which E(x) becomes the cut-off energy. As N(A,, x), 
m({x), E(x) and 7(£) are known, a(x) can be determined so as to satisfy Eq. (A, 29). 
The results obtained by such a procedure are shown in Fig. 1 and Table I, in case 
x=150 and 300. 

In case *«=1030, although FE, cannot be known directly from the observation, it 
seems to be reasonable to assume that the over-all multiplicity is similar to .that in 
the case where x is small, and that there also exists the critical energy E,. If we 
adopt as E, the largest value of E for which m(E, 1030) can be obtained from the 
observation and assume that m(E, 1030) is also given by Eq. (A, 28) in the energy 
region where E=E., a'(1030) can be determined so as to satisfy Eq: (A, 29) and gives the 
upper limit of a(1030). The over-all multiplicity m(E, 1030) obtained by such a 
procedure is shown in Fig. 1 and Table.I. On the other hand, there is another way 
to determine m(£, 1030) as follows; a(x) in Eq. (A, 28) increases with the increase of 
x, so that if @(300) is adopted as a(1030), this gives the lower limit of a(1030). In 


such a case, E,(1030) and (1030) can be determined so as to satisfy the following 
equation, 
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Ee(1030) 
Na 1030) ={ m (E, 1030) i(E) dE + | m(1030)| E 
a 


1 a(300) 
E,(1030 E,(1030) | *7(E)dE, (A, 30) 


where the over-all multiplicity obtained from the observation and its prolongation 
are used as m(E, 1030) in the 1st term on the right hand side. The over-all multi- 
plicity thus obtained is shown in Fig. 1 and Table I with the notation m™/ (E, 1030). 
The variations y’(A, 1030, 4E, 0) and y"(4, 1030, 4E,0) obtained by using m!(E, 1030) 
and m7(E, 1030) give the lower and upper limits of y, respectively. 
(2) In case of the neutron component 

The latitude dependence curves of the neutron intensity used to obtain the 
over-all multiplicity are those of the Simpson observation®@%? made at the atmos- 
pheric depthes x=305 and x=680. As these intensities are omnidirectional, the vertical 
intensities at x=305 can be obtained after the Gross transformation“) by using the 
absorption mean free path obtained by J.A. Simpson.” But, in case x= 680, it is not 
necessary to make the Gross transformation“®) and here the latitude dependence 
curve normalized to unity at equator is adopted. Assuming that the over-all multi- 
plicity is also expressed by Eq. (A, 28) in the high energy region, we can determine 
a(x) in the same way, and the results are shown in Fig. 2 and Table I. The over-all 
multiplicities obtained here are different in form with those obtained by S.B. Trei- 
man.“ This difference will be caused by the following facts that (1) the primary 
spectrums used are different from each other and (2) it is unreliable to determine m 
at the energy region in the neighbourhood of 10 Bev. by observations only. But the 
variation y at high latitude is almost not influenced by the adoptation of different 
over-all multiplicities, because the contribution of the high energy primary particles 
to the neutron intensity at high latitude is negligibly small as the latitude effect of 


the neutron component is very large. 
(Read Oct. 30, 1953) 
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Self-Reversal of Thermo-Remanent Magnetism 
of Igneous Rocks (III) *t 


By T. NAGATA, S. AKIMOTO and S. UYEDA 
(Geophysical Institute, Tokyo University) 


Abstract 


Besides the dacitic pumice of Mt. Haruna, an igneous rock 
having the tendency of self-reversal of thermo-remanent magnetiza- 
tion has been found in a pitchstone of Mt. Asio. Although the 
self-reversal property of the Asio ferromagnetic minerals is weaker 
than that of the Haruna ferromagnetic minerals, the magnetic and 
crystallographic properties of the former are fundamentally same as 
those of the latter. 
By thermo-magnetic separation, an ensemble of ferromagnetic 
grains in both the Haruna and the Asio rocks were divided into A, 
B and AB groups. The results of thermo-magnetic measurement, 
chemical analysis, and X-ray analysis of these three groups showed 
that the A-constituent is a titanomagnetite of the crystal structure of 
inverse spinel type having the Curie-point at 430°-510°C, and the 
B-constituent is an ilmenite-hematite solid-solution of the rhombo- 
hedral crystal structure having the Curie-point at about 230°C, while 
the AB grains are composed of both A- and the B-constituents. 
The measurement of thermo-remanent magnetization of a large 
number of single grain gave the result that only the AB grains can 
have the reverse thermo-remanent magnetization, the A- and the B- 
grains getting the ordinary normal thermo-remanent magnetization. 
§ 1. Tutroduction 
In the previous papers [1, 2], experimental evidence of the self-reversal of 
thermo-remanent magnetism of a dacitic rock of Mt. Haruna was reported: that is, 
the Haruna rock and the ferromagnetic minerals separated from it have the thermo- 
remanent magnetization, the direction of which is opposite to that of the magnetic 
field applied during their cooling from the Curie-point. It was also found out that 
the ferromagnetic minerals in the Haruna rock are composed of the group of ferro- 
magnetics similar to magnetite having the Curie-point at about 500°C and another 
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group of ferromagnetics having the Curie-point at about 230°C ; the former was called 
the A-constituent and the latter the B-constituent. The saturation magnetization of 
the B-constituent increases so linearly with decrease in temperature that this could 
not be identified to some of titanomagnetites which are usual sources of ferromagne- 
tism of rocks. 

Each grain of ferromagnetic minerals in that rock is composed of the A- and 
the B-constituents with their various ratios, the most grains being almost the A- 
constituent, while some of them are either the B-constituent or mixtures of A and B. 
It was suggested that only the grains consisting of both the A- and the B-constituents 
are responsible for the self-reversal of thermo-remanent magnetization; pure A- 
grains and B-grain can have the normal thermo-remanent magnetization. 

Thus it seems likely that the self-reversal of thermo-remanent magnetism of 
the Haruna rock is due to the mechanism proposed by Néel [3]; ¢.e. during a process 
of cooling in a weak magnetic field, the magnetic domains of the B-constituent are 
fixed in the direction of the strong demagnetizing field owing to the A-constituent 
magnetization which is parallel to the applied field, the intensities of both of the 
fixed A and B domains increasing with decrease in temperature ; the direction of the 
resultant magnetization of the sum of these two kinds of magnetization becomes 
opposite to that of the applied magnetic field below a certain temperature, since the 
rate of increase in the spontaneous magnetization with decrease in temperature is 
much larger for the B-constituent than for the A-constituent. Each of the above 
processes was experimentally verified for the case of the Haruna rock. We may 
therefore say that the Néel’s two-constituent hrontbone is valid for the reversed 
remanent magnetization of this rock. 

On the other hand, a large number of reversely magnetized igneous rock- 
bodies have been found hitherto at various localities over the world, and characteristics 
of thermo-remanent magnetism of some of these rocks were examined. For example, 
many samples of reversely magnetized rocks in Auvergne of France were examined 
by Roche, [4, 5], and his result shows that the thermo-remanent magnetization of 
these rocks is normal, ze. its direction is in agreement with that of the applied 
magnetic field. Similar studies were carried out by Hospers [6] for the basaltic 
lavas of Iceland which are reversely magnetized [7], and also by Vincenz [8] for the 
reversely magnetized dykes in England and Scotland. 

All of their results show that the thermo-remanent magnetism of these rocks 
is normal. We have also continued to examine the thermo-magnetic properties of 
reversely magnetized rocks in Japan, but we have been able to find out only two 
examples of igneous rocks having the self-reversal properties of thermo-remanent 
magnetism. 

The above facts would suggest that the typical self-reversal of thermo-remanent 
magnetism, which is reproducible in laboratory experiments, may be not common 
property of the reversely magnetized rocks. This was pointed out also by a number 
of scientists at the joint meeting of Royal Astronomical Society and the Geological 
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Society for the geophysical discussion on the magnetization of rocks [9]. 

Under such a circumstance, main problems in the study of the reversed thermo- 
remanent magnetism of rocks will be as follows: What parts of a large number of 
reversely magnetized rocks zz situ have the capability of the self-reversal of thermo- 
remanet magnetism: What are the exact conditions for the capability of the self- 
reversal in the physical, chemical and mineralogical properties of the ferromagnetic 
minerals in rocks. 

It will be quite reasonable to assume that various properties of many igneous 
rocks at present may have been appreciably changed from the corresponding original 
properties which these rocks had when they were produced in old geologic times. 
In other words, the thermo-magnetic properties of rocks at present might be appre- 
ciably different from those of the same rocks at their original fresh stage. As pointed 
out by Graham [10], the magnetic properties of rocks are not always _reproducible, 
even in the sense of rough approximation; that is, the magnetic properties of some 
rocks are so changed by heating only once up to a certain temperature that the repro- 
duction experiment of their thermo-remanent magnetism is hardly possible. This 
fact has long been noticed since it was found that the change in magnetization in 
some rocks is irreversible with respect to temperature [11, 12]. 

These facts may indicate that knowledge of the exact requirements for the 
capability of the self-reversal of thermo-remanent magnetism as well as the physical 
nature of general thermo-remanent magnetism itself will be able to lead us to the 
reliable conclusion about the cause of reversed magnetization of igneous rocks zz szfu. 

In this report, therefore, the result of experimental studies carried out as in 
detail as possible along the above-mentioned line will be summarized. 

§ 2. Properties of Thermo-Remanent Magnetism of Dacitic Pitchstone of Mt. Asio 

It has already been reported [1, 2] that, among a large number of the rock 
samples examined so far in our laboratory, the dacitic pumice of Mt. Haruna shows 
the typical capability of self-reversal of thermo-remanent magnetism. Similar studies 
have been extended in order to find other examples of the rocks having the same 
characteristic. 

What we have found so far is only one example of the dacitic pitchstone of 
Mt. Asio in Japan. The total thermo-remanent magnetization of this rock is by no 
means the reversed one, but a certain part of its partial thermo-remanent magnetism 
has the property of the self-reversal just similar to that of the Haruna rock. 

Fig. 1 illustrates the partial thermo-remanent magnetism (P.T.R.M.) of the 
ferromagnetic minerals separated out from the Asio dacitic pitchstone, where the 
temperature range 4T is kept 100°C. As clearly seen in this figure, P.T.R.M. 
produced by cooling from 300°C to 200°C in a magnetic field H, 72 », is reversed. 
The dependency of this J%},, upon H is illustrated in Fig. 2, where it is noticed that 
P.T.R.M. is reversed provided that H is less than the critical value 3.5 Oe., while it 
becomes normal if H exceeds this critical value. The general tendency of the curve 
Ji.1—H of the Asio rock is similar to that. of the Haruna rock described in the 


Self-Reversal of Thermo-Remanent Magnetism of Igneous Rocks (III) 171 


-4 


Fig. 2. Dependency of the partial thermo- 
remanent magnetization, J 39, 7 upon the 


Fig. 1. Partial thermo-remanent magnetism of external magnetic force, H. 
the ferromagnetic minerals in the Asio 


rock. (H=2.030e.) 


first report [1], except the critical value of H which is about 7.0 Oe. in the latter. 
In both the Haruna and the Asio rocks, the thermo-remanent magnetism of 
their ferromagnetic minerals is composed of the normal P.T.R.Ms. for the high 
(>400°C) and low (<200°C) temperatures and the reverse P.T.R.M. for the inter- 
mediate temperatures, the reversed component of the thermo-remanent magnetism 
taking place mainly during the temperature range from 300°C to 200°C. The 
remarkable discrepancy of the Asio sample from the Haruna sample is the magnitude 
of the reverse P.T.R.M.; ze. in the former, the reverse P.T.R.M. is so small in 
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Fig. 3. Specific susceptibility of the 
ferromagnetic minerals in 
the Asio rock. (H=1.350e.) 


comparison with the normal P.T.R.M., that the total 
thermo-remanent magnetism is normal. 

Thus, it seems likely that the fundamental 
properties of thermo-remanent magnetism of the 
ferromagnetic minerals of the Asio rock are mostly 
similar to those of the Haruna samples; the degree 
of the self-reversal is weakened in the former. 

Fig. 3 shows the change in; specific magnetic 
susceptibility with temperature of the Asio sample. 
This —T curve, forming the stepwise change, 
resembles that of the Haruna sample. In both 
cases, the —T curves have two Curie-points, ze. 
at about 540°C and about 280°C. This fact may 
also suggest that the ferromagnetic mineals of the 
Asio rock have the constitution nearly similar to 
that of the Haruna ferromagnetic minerals. 

§ 3. Thermo-Magnetic Separation of Ferromagne- 
tic Minerals 

It was decribed in the second report [2] that 
the preliminary trial to classify the ferromagnetic 
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minerals with respect to their Curie-points was fairly successful. The Peas 
procedure to separate ferromagnetic and non-ferromagnetic minerals at various 
temperatures will be called the thermo-magnetic separation. If the group of the 
ferromagnetic grains in a rock-sample is composed of various kinds of grains having 
different Curie-points, (9), we could find the spectrum /(@) of this group of grains 
with respect to their Curie-points with the aid of the thermo-magnetic separation. 

Let us assume that a grain is sufficiently ferromagnetic below 9 and completely 
non-magnetic above 9. Then the mass of ferromagnetic minerals which are non- 
magnetic at T; and ferromagnetic at JT; (4T=7;—T;>0), is given by 


Mj’ 7(6)48, (1) 


where M denotes the total mass of the whole ensemble of grains. If 4T=7;—T; is 
made sufficiently small, we can put that the mean value of f(9) for the range of 4T 
is nearly equal to f(®@) at 9=(T;+T7,)/2. 

In the practical experiment, all ferromagnetic grains at a temperature T are 
separated out, and then the temperature is decreased to T—4T, at which the same 
procedure of separating out the ferromagnetic grains is carried out. Repeating the 
same procedure at Successively lower temperatures, we .can get M\_ f(@d0,) 

Mi £949, mi " f(9) 49>. abe tant. 

Tf we take T= 600°C, i; f(9)d®@ is practically zero in most cases. This may 
indicate that the most igneous rocks contain magnetite and the other ferromagnetic 
minerals having the Curie-point less than that of magnetite. 

The principal parts of a thermo-magnetic separator are schematically shown in 
Fig. 4. As shown in this figure, the separation of ferromagnetic minerals is done in 
a N,-gas electric furnace in 
order to avoid the oxidation of 
samples. The grains to be 
separated are spread over in 
the form of single layer on a 
brass plate (B.P.); the laterally 
movable electro-magnet (E.M.) 
Sweeps over the plate and 
pulls up the ferromagnetic 
grains’ at T; the pulled up 
grains fall down into the trap 
by cutting off the electric 
current through the magnet; 
this process is repeated until 
no grain is pulled up by the 


Trap 


Fig. 4. Schematic view of thermo-magnetic separator. 
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Specimens is kept uniform within the error of 5°. However, the magnetic separation 
by means of the above method is subject to the uncertainty due to the intensity of 
magnetization and mass of the grains to be separated out. Let mass and the 
specific intensity of magnetization of a grain at temperature 7 under the magnetic 
force H of the magnet be denoted by m and J(H, T). Then, the necessary condition 
for pulling up the grain by the magnet is 


mH, T) mg, (2) 


where 0H/0z is the vertical gradient of magnetic force H on the plate on which the 
grain lies. If we put, 


¢/SH—c, (3) 


the critical condition for (2) is written as 


where J, denotes the critical temperature for the magnetic separation. Since H= 
1600 Oe. in the actual equipment. J(H, Tc) may be nearly equal to the saturation 
magnetization J,(T;) of the same sample at Ty. Hence (4) can be written as 


IAT c)=C. (4) 


The magnitude of C defined by (3) is an instrumental constant, and it may rule the 
_resolving power of the equipment for the thermo-magnetic spectrum f(9).. By using 
the small grains of pure nickel, the constant C was calibrated to be about 2.1 in the 
writers’ equipment. This means that the critical temperature T, of the magnetic 
separation is more or less lower than the Curie-point, since a grain is pulled up by 
the magnet only when its magnetization exceeds 2.1 e.m.u. per gram. 

Provided that dJ/(T)/dT near the Curie-point is sufficiently large, the difference 
49 between 9 and T; is approximately given by 


10 =6-Te=—C/( J), - (5) 


For the case of pure nickel grains 40 is about 4°. Since (0/,/0T)r.© of most ferro- 
magnetic minerals in rocks is of nearly the same order as that of nickel, 40 for 
them will be generally less than 5° degrees. It must be noticed however that To 
becomes appreciably different from 6 for such specimens as hematite whose J, and 
| 0J,/0T | are very small. 

Figs. 5 and 6 illustrate the thermo-magnetic spectra of the ferromagnetic 
minerals of the Haruna and the Asio rocks respectively. In these diagrams, the 
area of each rectangle represents ms,’ f(Tc)dT,=Mf(T,)4T, where 4T=T;—T; is 
not always constant in Fig. 6. In both spectra, there is a sharp increase of f(T;) at 
about T;~=450°C, and another small peak is seen at about 230°C. Hence we may 
conclude that the ferromagnetic minerals in the Haruna and the Asio rocks are 
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separation temperature of  ferro- separation temperature of  ferro- 
magnetic minerals of the Haruna magnetic minerals of the Asio rock. 

rock. Total mass M=1.015¢g. Total mass M=0.787 g. 


classified either to the group having Curie-point around 450°C or to the group of 
Curie-point of about 230°C. The former will be tentatively called the A-constituent 
group and the latter the B-constituent group. 

Strictly speaking, however, there may arise a question whether Ty really 
represents an approximate value of 9, if some non-magnetic substances are contained 
in the grain of ferromagnetics. 

Let the mass of non-magnetic substance at temperature T be denoted by mz,. 
Then the necessary condition for pulling up a grain composed of m, of non-magnetic 
substance and m, of magnetic one whose magnetization is J,(T) at T is given by 


mala TG — (mat) g=0 , (6) 


Putting m,/m,s=x, we get from (3), (4) and (6) that 
JAT,)=(1+)C. (7) 


This indicates that T, is lowered by x. If x is not large, 40 =0-—T, for this case is 
approximately expressed by ; 


40 =—(1+a)C/( HA) (8) 


T=8 


Since it is practically difficult to Separate exactly the ferromagnetic minerals 
alone from their mother rock, the above mentioned effect may probably take place 
in the results illustrated in Figs. 5 and 6. In other words, it may be ambiguous 
whether the spread of the spectral lines in those T, diagrams represents the real 


spread of the Curie-points or the spread of 40 owing to various different amounts 
of x for different grains. 
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§ 4. Change with Temperature of Magnetization of A, B, and AB-Grains 

The magnetization of each single ferromagnetic grain, which was separated out 
by the thermo-magnetic separation method, was measured in a fairly strong magnetic 
field at continual temperatures from 20°C to 600°C by a magnetic balance [13]. 

The magnetic field for measuring the magnetization was about 2400 Oe. Ac- 
cording to the measurement of the magnetization curve of these specimens, the 
intensity of magnetization observed in. H=2400 Oe. is not much less than that of the 
saturation magnetization. 
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Fig. 8. J(T)-T curve of single grain and 
thermo-magnetic separation spectrum 


: of the Asio rock. 
Fig. 7. J(T)-T curve of single grain and 


thermo-magnetic separation spectrum 
of the Haruna rock. 


The dependency of the nearly saturated magnetization (J(T)) upon temperature 
for single grains having various values of magnetic separation toeipetsnire (Tc) of 
the Haruna and the Asio rocks are illustrated in Figs. 7 and 8 respectively, where 
the abscissa represents the ratio of the intensity of magnetization J(T) of a single 
grain at T to that of the same grain at 20°C. Throughout results obtained, the 
curves J(T)-T are reversible with respect to temperature; ze. the J(T)-T curve in 
the cooling process is practically coincident with that in the heating process. 

As clearly seen in Fig. 7, the J(T)-T curves in one group have sain Curie- 
point at temperatures from 430°C to 510°C, taking the forms of the Curie ach 
The J(T)-T curves of another group have the Curie-point at about 230°C, their form 


176 TT. Naaara, S. Aximoro and S. Uvepa 


being almost a straight line representing a linear relation between J/(T) and T. 
Generally speaking, the former group of J(7)-T curves corresponds to the A-cons- 
tituent group defined in the preceding section, and the latter group to the B-constituent 
group. The specific intensity of magnetization of the A-constituent at 0°C is about 
80 e.m.u., while that of the B-constituent is about 20 e.m.u.. These values do not 
differ much from their respective saturation values. 

There are also several cases that the J/(7)-T curve seems to consist of the 
straight line in the temperature range below 230°C and the Curie curve for the 
temperatures above 230°C. There will be no room for doubt that the grains having 
the above properties of their J(T)-T curves consist of the A-constituent phase and 
the B-constituent phase. These grains which are composed of both the A- and B- 
phases will be called the AB-grains. The general tendency of the J/(T)-T curve of 
the Asio grains illustrated in Fig. 8 is almost the same as that of the Haruna 
samples ; z.e. these grains may also be devided into the three main classes: A, B, 
and AB. 

It seems likely that the magnetization J,, of the AB-grain in a strong magnetic 
field is expressed as 


Jan= Jat, (9) 


were m=m,+m,. Then, judging from the J(T)-T curves of the AB-grains, we may 
say that the ratio m,;/m, can take various values from zero to infinity in both the 
Haruna and the Asio rocks. For example, with the aid of (9), the ratio m2,/m, 
is estimated to be about 0.7, 31 and 235 for the three typical curves of the AB type 
illustrated in Fig. 7. 

However there is practically a limit for the value of the ratio which could be 
detected by measuring the J(T)-T curve. The detectable value of m,/m, in the 
authors’ equipment is in the range from about 0.5 to about 300. That is to Say, 
a very small content of the A-constituent within the AB-grain can rather easily be 
detected, but the B-constituent in an AB-grain can be detected only if its content 
exceeds a half of the whole grain. 

It will be concluded from the diagram in Fig. 7 that the Curie-point 9 of the 
A-constituent spreads over the temperature range from 430°C to 510°C and that of 
the B-constituent over the range from 230°C to 240°C. This spread of the Curie- 
point may indicate that neither the A- nor the B-constituent group is composed of 
ferromagnetic grains having a common constant structure, but each of them has 
a finite width spectrum with respect to the Curie-point. Let us now compare these 
groups of J(T)-T curves with the result of the thermo-magnetic separations of the 
same samples. The lower half in Figs. 7 and 8 represents the frequency distribution 
f(Tc) of the Haruna and the Asio grains respectively with respect to their tempera- 
ture TZ; of the thermo-magnetic Separation. Then it may be said that the spectrum 


t(Tc) can approximately represent the spectrum of Curie-point 4(0) except for the 
case of the AB- grains. 
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Since J,(T) is practically zero when T exceeds @,, the discussions summarized 
by (6), (7) and (8) can be applied to the case of the AB-grains provided that T>60,. 
Therefore, when m,;/m, in the AB-grains is larger, Tg for these grains will be the 
more lowered from 0,4. When Ty becomes less than @ x, the critical condition for the 
thermo-magnetic separation must be given by 


Tas To) = FAT (Te) +22] x To) = C (10) 


For example, the grain of 0,= 
470°C, @,=220°C and ms5/m,4=0.7, 
which is shown in Fig. 7, was separated 


out by the thermo-magnetic separator 
in the temperature range of 450°C- 
460°C. That the AB-grains are the 
aggregates of the A-constituent and 


the B-constituent was more directly 


proved as follows. The group of 
ferromagnetic minerals of 440°C<Tc< 


: i 4 8 i 

450°C was separated out by the thermo- B ———>  m/mem) A 
Fig. 9. Thermo-magnetic separation temperature 
(represented by a full thick line) with the 


rock; the grain-size is on average ratio of the content of A- and B-phases. 


magnetic separator from the Haruna 


about 0.3mm in diameter: then this 

group of grains was pulvurized into the small grains of 3~5 » in the mean diameter ; 
these small grains were again thermo-magnetically separated. The result showed 
that these small grains could be divided into the three groups: ze. 75.5% for 500°C 
<T,¢<525°C, 13.89% for 200°C< T<250°C, and the remaining 8.6% for non-magnetic 
minerals. The first and the second groups correspond to the A- and the B-constituents 
respectively. Thus, it will be clear that any ferromagnetic grain in the Haruna and the 
Asio rocks is composed of the A- and the B-constituents, which are distinctly separated 
out from each other in regard to their Curie-point. The above argument includes 
the special case that a grain consists of the A- or the B-constituent alone. 

It seems difficult to define uniquely the A- and the B-constituents, since the both 
constituents may have a finite range of their Curie-point. But the conception of the 
presence of the two constituents, A and B, in the ferromagnetic minerals of the 
Haruna and the Asio rocks will be acceptable. 

§ 5. Thermo-Remanent Magnetization of Single Grains 

The total thermo-remanent magnetization (TRM) of each single grain produced 
by cooling from 600°C in the geomagnetic field was measured by an astatic magneto- 
meter. The purpose of this experiment is to look for the source of the reverse 
thermo-remanent magnetization of the Haruna and the Asio rocks. This experiment 
was first done for 16 Haruna grains whose J(T)-T curves were obtained. 
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Each grain of the A-constituent having the Curie-point at 430°-510°C has weak 
normal TRM and also each single grain of the B-constituent having the Curie-point 
at 230°-240°C has weak normal TRM. The AB-grains which consist, as proved in their 
J(T)-T curves, of both the A- and the B-constituents have generally a strong reverse 
TRM. Some grains which have the Curie-point at 230°-250°C and so look like the 
B-constituent grains could have also a fairly strong reverse TRM. It seems likely 
therefore that these B-like-grains contain a small amount of the A-constituent which 
has effective influence upon the fixing of the magnetization of B-constituent, and so 
these will be called the aB grains. 

Thus it may be said that the ferromagnetic mineral which contains both the 
A- and the B-constituents within each single grain is mainly responsible for the 
reverse TRM. This result will indicate positive evidence of the two-constituent 
theory proposed by Néel [3] and a little modified by the authors [1]. 

In order to confirm the above argument, it will be desirable to carry out the 
similar examination of a Jarge number of ferromagnetic grains in the Haruna and 
the Asio rocks. Since the exact measurement of the J(T)-T curves of many grains 
are too laborious, it was assumed that the procedure of thermo-magnetic separation 
is able to indentify approximately the A-, the B- and the AB-grains. As illustrated 
in Figs. 5 and 6, the f(T¢) spectrum can be divided into two main classes; the group 
A at temperature higher than 400°C consists of the A-constituent and the AB-cons- 
tituent of m,/nz4<20, the latter being denoted here by Ab; the group B at tempera- 
ture lower than 230°C consists of the B-constituent and the AB-constituent of m,;/m4, 
>30, the latter being denoted by aB. The grains having JT; at temperaturs between 
230°C and 400°C will be denoted by ab, in which 30>,/m,>20, as estimated in the 
diagram shown in Fig. 9. 

There is no ab grain in the Haruna rock, as shown in Fig. 5. The grains of 
the group A have either the normal TRM of weak intensity or the reverse TRM of 
moderate or weak intensity; the first and the second cases may correspond practically 
to the A grain and the Ab grain respectively, where the AB grains containing rather 
small amount of the B-constituent will be practically considered to belong to the A- 
grains. The grains of the group B have either the normal TRM of weak intensity 
or the reverse TRM of strong intensity; the first may correspond to the B-grain and 
the latter to the aB grain. Among 89 grains of the group A, 67 are of the A type 
and 22 are of the Ab type; and among 17 grains of the group B there are 15 of the 
aB type and 2 of the B type. 

A similar examination was done for the Asio grains; 9 among 10 grains of the 
group A are of the A type, the remaining one having the type of Ab, while 3 of 5 
grains of the group B is of the aB type and the remaining two are of the B type. 

It may therefore be concluded that each single grain of the AB-constituent can 
generally have the reverse TRM, and that the thermo-remanent magnetization of 
a rock sample as a whole can be reverse when the sum of the reverse magnetization 


of these grains exceeds that of the normal magnetization of the other grains. 
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§ 6. Chemical Composition and Crystallographic Structure of A- and B-Grains 
Our problem will now be concerned with what minerals are the A- and the 
B-constituents in the Haruna and the Asio rocks. The result of chemical analysis of 
the A- and the B-constituents of the Haruna sample, made by courtesy of J. Ossaka, 
is given in Table I, and they are represented by the points A and B on the diagram 
of the ternary system of FeO-Fe.0,;-TiO, in Mol percentage shown in Fig. 10. In 
this diagram, the points Fe,0;, TiFeO; and Fe,O, represent respectively hematite, 
ilmenite, and magnetite, and the line connecting Fe,;O,; and TiFe,O, represents the 
complete solid solution between Fe,O, and TiFe,O,, being usually called the titano- 
TiO, magnetite. According to Pouillard [14], the solid 
solution of hematite and ilmenite is possible only 
in the two limited ranges of Fe,O;:TiFeO;, which 

are represented by two limited lines in Fig. 10. 
Then it will be obvious in Fig. 10 that the 
chemical composition of the A-constituent is nearly 
a titanomagnetite, which is not much different from 
the pure magnetite, and that the chemical composi- 


<- a) 0 tion of the B-constituent is nearly that of a solid 


ie solution of TiFeO, with Fe,O;; z.e. an ilmenite- 
Fig. 10. Chemical composition of A- 


and B-constituents in the 
Haruna rock, represented on In the last column of Table I, the chemical 
a FeO-Fe,0;-TiO, diagram 
in Mol. percentage. 


rich ilmenite-hematite solid-solution. 


composition of .the whole ferromagnetic minerals of 
the Haruna rock is given for comparison. Com- 


Table I. Chemical Composition of ‘‘A’’ and ‘‘B”’’ Constituents of the 
Ferromagnetic Minerals of the Haruna Rock. 


(By J. Ossaka and T. Katura) 


A B Total 


Weight % Mol % Weight % Mol % Weight % | Mol % 


Fe,03 58.7 39.4 43.5 26.9 57.9 38.7 
FeO 34.2 51.1 24.2 33.2 34.7 51.5 
TiO: 7.0 9.4 32.3 39.9 7.4 958 


paring this chemical conposition with that of the A- and the B-constituents, we may 
say that the largest parts of the Haruna ferromagnetic minerals are the A-constituent, 
the ratio of content of B to that of A being of the order of 1/100. This is in good 
agreement with the result of the thermo-magnetic separation dealt with in § 4. 

For the purpose of determining the crystal structure of these A- and B- 
constituents, their X-ray analysis was carried out with the aid of a “NORELCO” 
X-ray spectrometer. Examples of the X-ray spectrograms of the A-, the B- and the 
AB- grains of the Haruna rock are illustrated in Fig. 11. Here the grains of Tc>500°C, 
those of T7c¢<250°C, and those of 250°C < Tce<400°C were assumed to represent the A-, 


the B-, and the AB-grains respectively. 
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It was found that the A-constituent has the cubic crystal structure of inverse 
spinel type which is the same as that of magnetite, while the B-constituent has a 
rhombohedral crystal structure, ze. the same structure as that of hematite and 
ilmenite. ‘The crystal structure and the lattice parameters of the A-, the B-, and 
the AB- constituents of both the Haruna and the Asio rocks are summarized in 
Table II, where those of almost pure magnetite and ilmenite, observed with the same 
apparatus, are also given for comparison. 

It will be evident that the A-constituent chiefly consist of a magnetite-rich 
titanomagnetite, since the lattice constant of the first is nearly the same as that of 
the second. This is in good agreement with the results of both chemical analysis 


and the thermo-magnetic separation. 
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Fig. 11. Examples of X-ray spectrogram of A-, AB-, and B- constituents by ‘‘ Norelco”’ 
with Fe K radiation. 


Table Il. Crystal Structure and Lattice Parameters of A, B, and AB Constituents 
of the Haruna and the Asio Rocks. 


Source Grain Crystal Structure | Lattice Parameters 
A Cubic (Inverse spinel) a. =8.40340.002 A 
B Rhombohedral arn =5.4800.002 A 
Haruna Sea =55°08’ +01’ 
F 1 = 8.4 y 
AB Cubic + Rhombohedral jin = Sa7atP et 
Xen = 55°08! +05’ 
Cubic (Inverse spinel) a, =8.397+0.004 A 
Asio B Rhombohedral (ed =5.483+0.004 A 
Open =55°01' +03! 


AB | Cubic+Trace of Rhombohedral a, =8.405+0.003 A 
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Magnetite Cubic (Inverse spinel) a, =8.397+0.003 A 
Tiesenite Rh ie »=5.540-+0.002 A 
ombohedral On = 54°44’ 02! 


The crystal structure of the B-constituent is rhombohedral, but its lattice 
parameters are appreciably different from those of pure ilmenite; that is, they are 
shifted towards those of hematite. Taking the chemical composition into considera- 
tion, we may conclude therefore that the B-constituent is chiefly composed of a solid- 
solution of ilmenite with hematite. 

It will be seen in Fig. 11 and Table II that the X-ray spectrogram of the 
AB-grains is composed of both the cubic structure lines of the A-constituent and the 
rhombohedral structure lines of the B-constituent. This result indicates that the 
AB-grains are mixtures of the A- and the B-grains. 

It may be safely said therefore that all of thermo-magnetic, chemical, and 
crystallographic analyses of the structure of the A-, B-, and AB-constituents have 
shown the consistent result. 

§ 7. The Ferromagnetic Mineral of B-Constituent 

It has already been shown in the previous papers [1, 2] that the reverse TRM 
of the Haruna rock at the atmospheric temperature is caused by the marked increase 
of the reverse magnetization with decrease in temperature of the B-constituent, 
which overcomes the normal magnetization of the A-constituent. It has further been 
found in the present study that this B-constituent is a kind of the ilmenite-hematite 
solid-solution. 

On the other hand, it has been believed that both hematite and ilmenite are 
by no means the ferromagnetic minerals, and it was proved [15] that at least hematite 
has the antiferromagnetic arrangement of the spins within the crystal. Then we 
may presume that a particular solid-solution of these two non-ferromagnetic minerals, 
ilmenite and hematite, might be able to have such a strong ferromagnetic property 
as that of the B-constituent. Hence the magnetic properties of the B-constituent 
may deserve some detailed investigation. 

As already shown in Figs. 7 and 8, the Curie-point of the B-constituent is about 
230°C and its saturation magnetization at 0°C was determined to be about 20 emu/gr. 
This saturation magnetization of this B-constituent changes almost linearly with 
temperature between the Curie-point and 0°C. The measurement of the magnetization 
was extended down to the temperature of the boiling liquid air (about —180°C), the 
result showing that its linear relation with temperature holds almost perfectly 
between 230°C and —180°C. 

By courtesy of Néel and Pauthenet, the saturation magnetization (/,) of the B- 
constituent in the Haruna rock was measured by them for a wide range of tempera- 


ture from the Curie-point to the temperature of the liquid helium. Their result is 
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reproduced in Fig. 12. As clearly seen in this 
diagram, the saturation magnetization J, in- 
creases linearly with decrease in temperature 
down to about 70°K, becoming about 54 
e.m.u./gr. at 0°K. This intensity of J, at 0°K 
is more than a half that of magnetite at the 
same temperature. The above mentioned 
linear change of J, with temperature may 
suggest that the B-constituent is a kind of 


the ferrimagnetic material. raat 
On ene, tier hand, Tine. saat) eae Fig. 12. Saturation magnetization Jz of 
hematite has usually feeble ferromagnetic the B-constituent of the Haruna 
properties has often been reported [16], and rock, H=19,500 Oe. (After Néel 
: and Pauthenet) 
Néel [17] attempted with a successful result 
to explain the origin of this feeble ferromagnetism by assuming the presence of some 
deviation of the arrangement of spins from their perfect antiferromagnetic order. 
It might therefore be possible to assume that the fairly strong ferromagnetic property 
of the B-constituent is due to a kind of ferrimagnetism caused by a certain disturbance 
in the antiferromagnetic arrangement of spins in the ilmenite-hematite crystal. If 
so, it seems likely that this ferrimagnetism is due to the substitution 


2 Fet*++<Fett+Tit**+, 


where Fe*t*+ and Fe** have the magnetic moment of 5 and 4 in Bohr magneton 
number while Tit+*+ has no magnetic moment. 

In the case of the B-constituent, its chemical composition can be expressed 
approximately by Ti:FesO; and the intensity of its saturation magnetization is 
represented by about 5/3 Bohr magnetons per Mol. Then taking into consideration 
the condition that ilmenite and hematite may have the perfect antiferromagnetic 
arrangements such as (Ti,;+*+* Fe, **) (Ti: +*+* Fe: **) O,-- and (Fe***) (Fet**+) O;-~ 
respectively, we might be able to presume a suitable ferrimagnetic arrangement of 
Fe***, Fe** and Tit*** ions in the two sites of the crystal so that the arrangement 
fulfils the observed coditions for J, and the chemical composition. The simple nearest 
arrangement will be (Fe, ++ Fez +++ Ti, ++*+) (Fe,** Ti: ****) O3--, which gives Tis Fes 
O;, in chemical comsosition and 2 Bohr magnetons in spontaneous magnetization. 

The above argument may still be somewhat ambiguous, since the accuracy of 
both the chemical composition and the spontaneous magnetization of the substance 
may be not sufficient enough. It will be desirable therefore to study systematically 
the magnetic properties of the TiFeO,-Fe.O, series. 

§ 8. Conclusion and Unsolved Problems 

A phenomenological mechanism of self-revercal of remanent magnetization of 
the AB-grains has already been proposed in the previous paper and summarized in 
§ 1 of this article. In the present study, it has been confirmed that the A-constituent 
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is a titanomagnetite having the Curie-point at 430°C-510°C and the B-constituent is 
a ferrimagnetic solid solution of ilmenite and hematite, having a rhombohedral 
crystal structure and the Curie-point of about 230°C. The observed large value of 
the temperature coefficient of Je below ©, in comparison with that of J, seems to 
be sufficiently large to verify the two-constituent theory of reverse TRM. 

However there still remains problem concerning detailed mechanism of the 
assumed magnetic coupling between the A- and the B-phases within each AB-grain. 
This is by no means self-evident. According to the result of microscopic observation 
of the surface of the AB-grains, the B-phase of ilmenite type of the form of a thin 
band sometimes intergrows within the A-phase. It seems that this behaviour of the B- 
phase is favourable to the assumption of strong anti-parallel coupling between Ja 
and J;. From the geophysical stand point of view, it will be desirable to know how 
often the necessary condition for the anti-parallel coupling between J, and J, can 
takes place within the actual ferromagnetic minerals. In other words, the possibility 
of the said coupling will be a significant problem in relation to the generality of the 
phenomenon of self-reversal of thermo-remanent magnetization. It seems therefore 
that this problem well deserves the further detailed study. 

Another important problem will be the degree of petrological commonness of 
the ilminite-like B-constituent in rocks. As clearly shown in this study, the amount 
of the B-costituent is only about one per cent of that of the A-constituent even in 
the Haruna rock which shows the remarkable reverse thermo-remanent magnetization. 
The presence of this particular ferrimagnetic phase in this rock had been overlooked 
even in the authors’ own study of the general magnetic properties of rocks until its 
characteristics of themo-remanent magnetism were revealed. 

On the other hand, the crystals of ilmenite, so-called in petrology, appear in 
almost all igneous rocks and further they are usually associated with magnetite 
crystals, as observed under a microscope. But very little work: has been done with 
respect to the detailed magnetic properties of ilmenite in rocks. It might be possible 
therefore to assume that some ilmenites have the ferrimagnetic properties such as 
those of the B-constituent. If the remanent magnetization of this ilmenite-like 
ferrimagnetic phase is stabler than that of the ordinary titanomagnetite, we may 
expect that the former magnetization alone has survived at present after a long life 
of rocks since they were formed in geologic time. If further there was the strong 
anti-parallel coupling between the ilmenit-like ferrimagnetics and the titano- 
magnetites, we may find that the rocks show the reverse magnetization im sztu, but 
we may not be able to reproduce the reverse magnetization by a simple laboratory 
experiment of thermo-remanent magnetism, except for the special case that J, of the 
ilmenite-like ferrimagnetics is stronger than J, of the titanomagnetites at the 
initial stage of the creation of thermo-remanent magnetization. 

The above-mentioned is merely an example of hypothetical possibility, just 
for emphasizing a significance of presence of the ilmenite-like ferrimagnetics in 
rocks. This problem must be examined in detail in future, 
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